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Abstract 

We present a systematic, and fully analytic, study of the ratio TZ p a between the 
gluon distribution in a nucleus and that in a proton scaled up by the atomic number 
A. We consider initial conditions of the McLerran-Venugopalan type, and quantum 
evolution in the Color Glass Condensate, with both fixed and running coupling. We 
perform an analytic study of the Cronin effect in the initial conditions and point out 
an interesting difference between saturating effects and twist effects in the nuclear 
gluon distribution. We show that the distribution of the gluons which make up the 
condensate in the initial conditions is localized at low momenta, but this particular 
feature does not survive after the quantum evolution. We demonstrate that the 
rapid suppression of the ratio 1Z p a in the early stages of the evolution is due to the 
DGLAP-like evolution of the proton, whose gluon distribution grows much faster 
than that in the nucleus because of the large separation between the respective 
saturation momenta. The flattening of the Cronin peak, on the other hand, is due 
to the evolution of the nucleus. We show that the running coupling effects slow down 
the evolution, but eventually lead to a stronger suppression in 1Z p a at sufficiently 
large energies. 
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1 Introduction 



Motivated especially by the recent experimental results for gold-gold (Au-Au) [1,2] and 
deuteron-gold (d-Au) collisions at RHIC [3-6], there is currently a large interest in the 
physics of collective phenomena in the wavefunction of an energetic nucleus, which could 
explain, for instance, why the gluon distribution of a large nucleus (with atomic number 
A ^> 1) is not simply the incoherent sum of the gluon distributions of the A constituent 
nucleons. For sufficiently high energies (and, marginally at least, also for the energies at 
RHIC), one expects these collective phenomena to be associated with non-linear effects 
favored by the high gluon density, which rises rapidly with the energy, and also with A 
[7-12]. Thus, these effects are expected to be more pronounced in a (large) nucleus as 
compared to a proton, and this difference should be particularly important at not so high 
energies, where the proton is still in a linear regime. 

This difference could explain some of the remarkable regularities observed in the results 
for d-Au collisions at RHIC, like the Cronin enhancement in particle production at cen- 
tral rapidity (y = 0)[3], and the 'high-/cj_ suppression' in the particle yield at 'forward 
rapidities' (y = 2 — 3 in the deuteron fragmentation region) [4-6]. Specifically, the Cronin 
peak [13] can be understood as the result of Glauber-like multiple scattering off the gluon 
distribution produced by uncorrelated "valence quarks" [14] (as described, e.g., by the 
McLerran-Venugopalan model [10], see also [15]), whereas the high-/cj_ suppression can 
result from non-linear effects in the quantum evolution with increasing energy [16]. (See 
also [17] for alternative interpretations of these results.) On the other hand it is most 
likely that a similar suppression seen in Au-Au collisions at central rapidity [2] is due to 
jet quenching through final state interactions [18]. 

More precisely, the experimental results are used to construct the ratio RdAu(k±,y) be- 
tween the number of particles produced in a d-Au collision (with a given transverse 
momentum k± and at a given rapidity y) and the corresponding number for a proton- 
proton collision scaled up by the number of collisions. The deviation of this ratio from one 
- like a Cronin peak (RdAu > 1 at intermediate momenta), or the high-Zc^ suppression 
(RdAu < 1 at generic, hard, momenta) — can be attributed to the difference between 
the gluon spectrum (or "unintegrated gluon distribution") in the nucleus, and that in the 
proton scaled up by A. In turn, this difference is measured by the ratio: 



n pA (k„y) = A1/3(pp{k± ^ yy (i.i) 

between the gluon occupation factors (see Eq. (2.1)) in the nucleus and the proton 1 . 



1 Our conventions are such that the more standard "unintegrated gluon distribution" is equal 
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Previous studies in the literature show that it is indeed possible to compute the cross- 
section for particle production in terms of the gluon distributions in the target and the 
projectile [7,14,19-22], and, moreover, any qualitative trend seen in the ratio (1.1) of 
the gluon distributions gets transmitted to the corresponding ratio for the cross-sections 
[16,23-25], which is the quantity relevant for RHIC. 

Motivated by this observation, and also by the conceptual importance of the ratio (1.1) 
as a direct measure of collective effects in the nuclear wavefunction, we shall devote this 
work to a systematic analysis of this quantity within the effective theory for the Color 
Glass Condensate (CGC) [26], which is the appropriate framework to describe non-linear 
effects in the gluon distribution within QCD. 

Although inspired by the RHIC phenomenology, our analysis does not aim to the descrip- 
tion of the data (for such attempts, see e.g. [27,28]), but is rather oriented towards the 
fundamental understanding of phenomena like the emergence of the Cronin peak in the 
ratio TZpA and the suppression of this peak with increasing energy. Our main objective is 
to present a calculation for such phenomena from first principles, and which is fully analytic. 
These requirements entail strong constraints on the kinematical range which is accessi- 
ble to our analysis, and also a significant loss of accuracy, due on the one hand to the 
limited accuracy of the perturbative framework that we shall use, and on the other hand 
to the additional approximations that we shall perform within this framework, to allow 
for analytic calculations. In particular, our focus on the gluon distribution, rather than 
on the cross-section for multiparticle production, is also motivated by the requirement 
of (analytic) calculability, together with our concern to avoid theoretical uncertainties as 
much as possible. 

As already mentioned, our theoretical framework will be the CGC effective theory[10, 26, 
29-36] (see also [11, 12] for recent review papers), which is a lowest-order (in a s ) formalism 
in QCD, although obtained after elaborate resummations: These resummations include 
only those diagrams of perturbative QCD in which the powers of a s are enhanced by either 
large logarithms of the energy (i.e., by powers of y ~ In s), or by high density effects (e.g., 
by the gluon occupation factor, which at saturation is of order l/a s ). Still, it should be 
emphasized that this formalism becomes more and more accurate with increasing energy 
because, first, it does include the leading effects at large y, and, second, the saturation 
momentum Q s {y) — the typical momentum of the gluons which make the condensate - 
increases very fast with y [7, 37-42] , so the relevant coupling constant decreases at high 
energy: a s (Q s (y)) oc 1/y^ at large y. Thus, the results that we shall obtain here represent 
the actual prediction of QCD for sufficiently large y. 

This discussion also shows that the inclusion of running coupling effects in the formalism is 
crucial in order for the asymptotic freedom of QCD to become operative at high energies. 

to SAfA(k±, y), with Sa the area of the nuclear disk. 



4 



The running of the coupling is a higher order effect in a s , so its inclusion in a leading 
order formalism is by definition ambiguous. But physical intuition about the typical scales 
involved in the interactions, and also the experience with BFKL equation [43] — which in 
the context of the CGC formalism is the linear equation describing the approach towards 
saturation, and is presently known to next-to-leading order in a s [44, 45] — , will permit 
us to effectively use a one-loop running coupling, with the scale for running set either by 
the gluon transverse momentum k±, or by the saturation momentum. 

As mentioned above, another source of accuracy loss are the further approximations 
needed to make the CGC formalism tractable via analytic calculations. Whereas the 
initial condition that we shall use at low energy — namely, the McLerran-Venugopalan 
model [10] - is sufficiently simple to allow for exact (analytic) calculations, the subse- 
quent evolution with increasing y is described by complicated non-linear equations which 
couple n-point correlations with any number of points n [26,32,34]. This is similar to 
the infinite hierarchy of Schwinger-Dyson equations, and within the CGC theory it can 
be summarized into a closed functional equation for the generating functional of the cor- 
relations [26,31,34], also known as the JIMWLK equation. Exact numerical solutions to 
this functional equation have recently become available [46], but these are still difficult 
to handle, so most numerical analyses (including that of the ^p^-ratio in Ref. [25]) have 
instead focused on a simpler non-linear equation, the Kovchegov equation [33], which 
is a kind of mean field approximation to the general JIMWLK equation [47,48], and is 
expected to work reasonably well for a large nucleus at not so high energies (for numerical 
analyses of the Kovchegov equation, see [49-52]). 

But whereas the numerical calculations are undoubtedly important, they cannot fully 
supplant the analytic studies as far as our physical insight and fundamental understanding 
are concerned. The numerical results in Ref. [25] have revealed very interesting features 
- notably, an extremely fast suppression of the ratio TZ p A{k±,y) at generic momenta 
with increasing y, and also the complete disappearance of the Cronin peak after only a 
short evolution — , but these observations have also raised new questions, since the precise 
mechanism behind such a rapid evolution has remained unclear. In particular, it was not 
clear whether the two effects alluded to above (the high-/cj_ suppression and the flattening 
of the Cronin peak) were caused by the same mechanism, or not. Also, the specific roles 
played by the proton and by the nucleus during evolution have not been elucidated. 

The previous analytic investigations in the literature appeared before the numerical results 
in Ref. [25], and cannot fully explain the latter. For instance, while the arguments in 
Refs. [16, 24] do explain the suppression at relatively large y and for high k±, they cannot 
explain why this suppression is so rapid in the early stages of the evolution, or why the 
Cronin peak flattens out so fast. Indeed, the approximations used in these previous studies 
do not apply for transverse momenta in, or near, the nuclear saturation region, nor for 
very small values of y. Besides, running coupling effects have never been considered in 



5 



the previous studies, either analytic or numerical. As we shall see, these effects lead to 
important changes, especially for the evolution at late stages. More generally, what seems 
to be still missing is a coherent, qualitative and quantitative, picture explaining what are 
the specific features of the MV model which imply the emergence of a well pronounced 
peak at low energies, and also what are the generic features of the quantum evolution 
which lead to the rapid flattening of this peak and to a general suppression in the ratio 
TZpA at all non-asymptotic momenta. 

Our subsequent analysis is intended to fill this gap, through a systematic study of the ratio 
7Z p a in the MV model and of its evolution with y, for both fixed and running coupling. 
Being systematic, our analysis has necessarily some overlap with previous studies in the 
literature, especially with Ref. [24], with which we share some results and conclusions, 
notably about the high-/cj_ suppression. In our opinion, what distinguishes the present 
analysis from such previous investigations, in addition to its aim to completeness and 
the treatment of the running coupling case, is the constant effort towards elucidating the 
fundamental reasons for the observed behavior, and also the use of a coherent scheme of 
approximations, which will be carefully justified, and whose limitations will be discussed. 

In what follows, we shall briefly describe the picture which emerges from these calculations 
when increasing y from y — 0, and at the same time present the structure of the paper : 

i) Initial conditions at y = : Cronin peak in the MV model (cf. Sect. 2) 

Although a Cronin peak in the ratio TZ p A(k±) has been clearly seen in numerical simu- 
lations of the MV model [22-25,53,54], its analytic study has been hindered so far by 
the non-linear aspects of the problem, which must be treated exactly. In Sect. 2, we shall 
present a complete, analytic, study of the nuclear gluon distribution in this model, which 
will enable us to compute the magnitude and the location of the peak, and clarify the 
conditions for its existence: The emergence of a well pronounced peak in the MV model 
reflects the peculiar redistribution of gluons in transverse phase space, under the influence of 
the non-linear effects responsible for saturation. In turn, this issue has two aspects: 

The very existence of a peak can be anticipated on the basis of very general arguments (like 
the sum rule introduced in Ref. [24]; see also Refs. [23,25] and Sect. 2.4 below), which 
reflect two basic properties of the MV model: (a) the fact that the gluon distribution 
saturates, at a value of order l/a s , for low momenta k± < Q S (A), and (b) the fact that, 
at large momenta k± ^> Q S (A), the nuclear gluon distribution is simply the incoherent sum 
of the individual distributions of the A nucleons (so that 1Z p A{k±) — > 1 when k± — > oo). 
As we shall see in Sect. 2.3, these two properties alone imply the existence of a peak at 
intermediate momenta. The height of the peak is of order pA = In Q 2 S (A) / Aq CD ~ In A 1 / 3 , 
and thus is parametrically enhanced at large A. 

Moreover, the fact that the peak is so pronounced, and located in the vicinity of the saturation 
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momentum Q S (A), is the consequence of the fact that the distribution of the saturated 
gluons which make up the condensate is compact, in the sense that it vanishes exponentially 
at momenta above the saturation scale. This feature has not been recognized in previous 
studies of the MV model. Of course, on top of this compact distribution there is also a 
power-law tail, which represents the sum of the 'twist' contributions usually mentioned 
when discussing the high-A;^ behavior of the distribution. But as we shall see, for momenta 
k± ^ Q S (A), the sum of the twist terms is parametrically suppressed at large A as compared 
to the compact distribution which represents the condensate. Because of that, the overall 
distribution has a rapid, exponential, fall-off at momenta just above Q S (A), leading to 
the pronounced peak seen in numerical simulations [22,23,25,53,54]. 

ii) y > : Non-linear quantum evolution in the CGC (cf. Sect. 3) 

In Sect. 3 we shall present our approximations to the solution to the non-linear equations 
which describe the evolution of the gluon occupation factor with y. 

In the saturation region at low momenta, physics is fully non-linear, but for k± <C Q S (A, y) 
the solution can be constructed in a mean field approximation [35,36]. The result is an 
universal, and slowly varying, function of z = k\/Q 2 s {A, y) : if a oc (l/a s ) \nl/z for z <C 1 
[35-37,55]. The solution in the transition region around Q s (A,y) is not known in general 
(see however Sect. 5.3), but we do know that, along the saturation line k± = Q s (A,y), tp^ 
is constant and of order l/a s . 

At high momenta k± ^> Q s (A,y), the gluon density is low, and the general evolution 
equations reduce to the BFKL equation [43], which, although linear, is still sensitive to 
saturation effects, through its boundary conditions [7, 37-41]. The solution to this equation 
in the presence of saturation has attracted much attention in the recent literature [38- 
41,48], where analytic approximations have been constructed for both fixed and running 
coupling. Here, we shall simply summarize the relevant results, and discuss their range of 
applicability. Still, a few steps in the construction of these approximations will be outlined, 
to clarify some subtle points and reveal the generic features of the evolution which will 
be needed in the subsequent discussion of the ratio 1Z p a- 

We shall quite generally distinguish between two physical regimes, which are described by 
different approximations: (I) a "DLA" ("double-log accuracy") regime at very large mo- 
menta k± ^> Q g (A, y), in which the evolution is dominated by the transverse phase-space 
p(A,k±) = \nk\/Ql(A), and (II) a "BFKL regime" within the range Q s (A,y) < k± < 
Q g (A,y), where physics is linear but influenced by saturation, and the gluon distribution 
preserves the geometric scaling [56] property characteristic of saturation: it depends upon 
the kinematical variables kj_ and y only through the ratio z = k^/Q^A, y) [38, 39]. Here, 
Q g (A,y) is the "geometric scaling momentum", which grows faster than the saturation 
momentum with y, and marks the upper bound of the geometric scaling region [38] (see 
also [39-41, 48]). With a running coupling, the discussion of the physical regimes becomes 
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slightly more involved, and will be presented in Sect. 3.3. 

The approximations described in this section will be sufficient to study the suppression 
in the ratio 1Z p A(k±,y) with increasing y at generic momenta. On the other hand, they 
are not sufficient to describe the flattening of the Cronin peak (see below). 

iii) y > : The general argument for W\gh-k± suppression (cf. Sect. 4) 

In Sect. 4 we shall identify and characterize the general features of the quantum evolution 
which are responsible for the suppression of the ratio 7Z p a at generic momenta. 

The rapid suppression of the ratio lZ p A(k±,y) with increasing y is shown to be the con- 
sequence of the strong dissymmetry between the quantum evolution of the nucleus and 
that of the proton, which in turn reflects the large separation of scales between the re- 
spective saturation momenta. Physically, it is so because, for the same values of y and k±, 
the transverse phase-space available for the evolution, namely p(A,k±) = \nk\/Q 2 s (A), 
is larger for the proton than for the nucleus (since Q S (A) ^> Q s (p)). In more technical 
terms, the proton evolves faster because of general properties of the kernel of the BFKL 
equation, like the convexity of its eigenvalue x(t)> which accelerate the evolution with 
increasing p. 

Moreover, the suppression rate dlnTZp^/dy is largest at small y and for not so large trans- 
verse momenta (as compared to the nuclear saturation momentum), since in this regime 
the difference between the evolution of the proton and that of the nucleus is most pro- 
nounced. This explains the rapid suppression observed in the early stages of the evolution 
in the numerical results in Ref. [25]: This is the consequence of the DGLAP [57] evolution 
of the proton (as described here by DLA), while the nucleus evolves comparatively little 
(except at extremely large momenta). 

The general properties of the evolution equations can also be used to show that, for fixed 
y with a s y ^ 1, the ratio TZ p A(k±,y) is monotonously increasing with k±. Thus, the Cronin 
peak flattens out during the first l/a s units of rapidity. 

Another interesting consequence of the evolution is a reversal in the A-dependence of the 
ratio 1Z p a when increasing y (see also Refs. [23,24]) : Whereas at y — 0, the height of the 
Cronin peak — and, more generally, the ratio lZ p A(k±) at generic momenta around Q S {A) 
- is (logarithmically) increasing with A, this tendency is rapidly reversed by the evolution: 
after only a small rapidity increase a s y ~ 7Z p A{k±, y) becomes a decreasing function 

of A for any k±. This property could be related to a similar change of behavior in the 
centrality dependence of the ratio RdAu measured at RHIC [4-6]. 

iv) The flattening of the Cronin peak (cf. Sect. 5) 
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The evolution of the proton alone produces a rather uniform suppression in the TZpA-ratio 
at all momenta. Thus, by itself, this cannot wash out any structure present in the initial 
conditions, like the Cronin peak. Therefore, the disappearance of the peak is necessarily 
related to the evolution of the nucleus. 

In order to study this phenomenon, we shall use the Kovchegov equation to compute the 
evolution of the nuclear gluon distribution in the first step Ay in rapidity, with a s Ay <C 1. 
This calculation allows us to follow the evolution of the peak at least up to the rapidity yo, 
with a s y ~ (In 2 pa)/pa 1, where the height of the peak becomes of order one. As the 
calculation shows, the effect of the evolution is to generate power law tails which progressively 
replace the original exponential tail of the gluon distribution at saturation. Because of that, 
the peak flattens out, and moves up to larger momenta. 

The flattening is related to the difference in the nuclear evolution at momenta below and, 
respectively, above the saturation momentum Q s (A,y). For k± < Q s (A,y), the gluons 
in the nucleus are saturated, and evolve only slowly. For k± > Q s (A,y), the nucleus is 
in the linear, BFKL, regime, and the corresponding gluon distribution increases much 
faster (although not so fast as the proton distribution at the same k±). Because of this 
dissymmetry, the Cronin peak gets tilted up, and flattens out very fast. Although we 
cannot control analytically this evolution until the complete disappearance of the peak, 
we shall check that the peak has flattened out when a s y ~ 1. 

v) High— suppression: the detailed picture (cf. Sects. 5 and 6) 

In Sect. 6, we shall give the detailed picture of the phenomenon of "high-A;_i_ suppression" 
in the kinematic plane y — In k\. (Previously, a similar analysis will be presented in Sect. 
5, but only for the evolution along the nuclear saturation line k± = Q s (A,y).) In fact, 
the suppression will be seen to occur at all momenta, and not only at those which are 
"high" in the sense of the present analysis (and which are such that k± ^> Q s (A,y)). 
Still, the suppression is stronger for momenta above the proton saturation scale Q s (p,y) 
- since this is where the proton evolves faster — but below the nuclear geometric scale 
Q g (A, y), since for even larger momenta, both the proton and the nucleus are in the DLA 
regime, and their respective evolutions almost compensate in the ratio 7Z P A(k±,y). Thus, 
not surprisingly, the ratio is found to asymptotically approach one from the below when 
increasing k± at fixed y. But since Q g (A,y) increases rapidly with y, it is clear that, for 
a s y ^ 1, a strong suppression will be visible at all non-asymptotic momenta. 

The explicit analysis in Sect. 6 will confirm the general trend of the evolution anticipated 
in Sect. 4: For generic momenta k± (such that the nucleus is either at saturation, or in the 
BFKL regime), the suppression in the ratio TZ p A{k±,y) is very fast in the early stages of 
the evolution, when the proton is in the DLA regime, but with increasing y it slows down, 
and eventually 1Z p a stabilizes at a very small value, proportional to an inverse power of 
A 1 / 3 . Besides, for fixed y (with a s y <^ 1), the ratio is monotonously increasing with k±, 
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from a very small value ~ l/A 1 ^ for k± ~ A QCD to 1Z pA ~ 1 when k± — > oo. 



y) Qs(A, y) 




Fig. 1. Physical regimes for evolution in the kinematic plane y — lnk 2 ^ (for fixed coupling, 
for dehniteness; see Fig. 10 for the corresponding picture with a running coupling). Both the 
saturation momentum Q S (A) and the 'geometric scale' momentum Q 9 (y) rise exponentially 
with y (for both the proton and the nucleus), and thus are represented by straight lines in 
the logarithmic scale of the plot. As visible on this plot, the geometric scale rises faster (its 
logarithmic slope is roughly twice as large as that of the saturation momentum). 

In order to describe this evolution in more detail, it is convenient to follow a line which 
is parallel to the saturation line in the kinematic plane y — In k\ (see Fig. 1). That is, we 
simultaneously increase y and k± in such a way to keep the ratio z = k"j_/Ql(A,y) fixed. 
If z = 0(1), then when increasing y from zero, the proton starts in the DLA regime, 
because Q g (p,y) <C Q s (A,y) at sufficiently small y. But with increasing y, Q g {p,y) grows 
faster than Q s (A,y), so at some 'critical' value y c the proton changes from the DLA to 
the BFKL regime. The value of y c is however different for fixed or running coupling, and 
so is also the evolution for larger rapidities y > y c , as we explain now: 

(a) With fixed coupling, a s y c ~ pa, and most of the suppression is achieved already for 
y ^ y c . For y > y c , the suppression slows down significantly because the BFKL evolutions 
of the proton and of the nucleus almost cancel in the ratio (1.1) (due to the fact that the 
respective saturation momenta evolve in the same way with y [42]). Thus, at large y, the 
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ratio TZ p a{z, y) decreases only slowly, and eventually stabilizes, when a s y ?t p\, at a small 
value TZpA ~ 1/ {A 1 ^pa) 1 ~ 1 ■, which shows a weak dependence on z. Here, 7 ~ 0.63 is the 
BFKL saddle point (or 'anomalous dimension') near saturation [38,39]. 



(b) With running coupling, y c is parametrically larger, y c oc p 2 A , but the suppression is 
pursued within a significant range of y above y c , because in that range the proton and 
nuclear saturation momenta evolve differently with y [42] : The proton saturation mo- 
mentum Q s {p,y) grows initially faster, and catches up with that of the nucleus when 
y ~ p\. For larger y, Q s (A,y) ps Q s (p,y), and the ratio stabilizes at the constant value 
lZ pA ~ I/A 1 / 3 which is even smaller than the corresponding limit for fixed coupling. In 
fact, this limiting value is simply the factor of A" 1 / 3 introduced by hand in Eq. (1.1). 
This reflects the universality of the gluon distribution which, within a wide range of k±, 
depends upon the hadron species only via the saturation momentum. 



It should be also stressed that, although the suppression in 1Z p a is eventually stronger 
with running coupling, the evolution leading to this suppression is slower in that case, in 
the sense that it takes a considerably larger interval of rapidity before the final limiting 
value is reached. Besides, all the intermediate stages of the evolution (like changing from 
the DLA to the BFKL regime with increasing y at fixed z) require a longer rapidity 
evolution with a running coupling than with a fixed one. 



To conclude this introductive discussion, let us emphasize that all the phenomena dis- 
cussed so far — the Cronin peak in the initial condition, its disappearance with increas- 
ing y, and the high-Zc^ suppression — are hallmarks of saturation, which however reflect 
different aspects of saturation, and of the quantum evolution towards saturation : The Cronin 
peak reflects classical saturation, i.e., the saturation via non-linearities in the classical field 
equations, whereas the color sources in these equations (the 'valence quarks') are uncorre- 
lated. (Alternatively, this reflects incoherent, Glauber-like, multiple scattering.) This is, of 
course, the content of the MV model, but it could also be a reasonable approximation for a 
large nucleus at not so high energies, where the effects of quantum evolution are negligible. 
Furthermore, the driving force towards both high-Zc^ suppression and the flattening of 
the Cronin peak is the linear evolution (either DGLAP, or BFKL). But notwithstanding 
this, these phenomena are still a signal of saturation, since they occur only because of 
the mismatch between various regimes of evolution, which is made possible by saturation : 
For the flattening of the peak, this is the mismatch between the nuclear evolution below 
and above the saturation scale, whereas for the high-Zc_i_ suppression, this is the mismatch 
between the (linear) evolutions of the proton and of the nucleus, which originates in the 
difference between the respective saturation momenta. 
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2 Cronin effect in the initial conditions 



2.1 Generalities 



By the "unintegrated gluon distribution" we shall more precisely understand in what 
follows the gluon occupation factor, i.e., the number of gluons of given spin and color per 
unit phase-space in a nucleus with atomic number A: 



( 27 3 ^N A 

^^ )E ^ w^ - (2 - 1} 

Here, y denotes the rapidity (related to the longitudinal momentum fraction x of the 
gluons via y = ln(l/x)), k± is the transverse momentum, and b± is the gluon position 
in transverse space. For simplicity, we shall consider a hadron which is homogeneous 
in the transverse plane, and we shall suppress the fe^-dependence in all the formulae. 
Within light-cone quantized QCD, the gluon occupation factor (2.1) is related to a gauge- 
invariant 2-point function of the color fields, which can be computed within the CGC 
effective theory [12,26]. More precisely, the evolution of this quantity with y can be 
computed within perturbative QCD (to the accuracy of the effective theory), but the 
initial conditions at y = remain non-perturbative, and require a model. It turns out 
that the choice of this model has a strong influence on the physical problems that we 
would like to address, especially as far as the existence of a Cronin peak [13] is concerned. 

Physically, the initial conditions that we are interested in correspond to a large nucleus 
(A ^> 1) in a regime of intermediate energies. The initial energy should be high enough for 
the gluons with the smallest values of x to be coherent with each other, but low enough for 
the quantum effects in the gluon distribution to be negligible. These conditions require 
x ^ 1/A 1 / 3 , but at the same time o s lnl/a;o <C 1, where x denotes the longitudinal 
fraction of the slowest gluons within the nucleus at the initial energy. For instance, for a 
gold nucleus and with a s ~ 0.2, these requirements can be satisfied by choosing xo ~ 1CT 1 . 
(From now on, the rapidity variable will be understood to represent the difference from 
this original rapidity: y = ln(xo/x).) Under these conditions, the gluon distribution in 
the nucleus is simply the result of classical radiation from the valence quarks. But if A is 
large enough — which is what we assume here — , the resulting gluon density can be still 
very high, and thus favor non-linear effects which lead to gluon saturation at sufficiently 
low k±, below some characteristic scale Q S (A). 

A simple model which encompasses these physical conditions is the McLerran-Venugopalan 
(MV) model [10, 12], in which the total color charge in the nucleus is the incoherent sum of 
the individual color charges of the valence quarks, and the corresponding saturation mo- 
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mentum scales like Q S (A) ~ A 1 / 6 (since Q 2 S (A) is proportional to the color charge squared 
per unit transverse area). If A is large enough, this scale is hard (Q 2 (A) 3> Aq CD ), and 
the matter made of the gluons (the CGC) is weakly coupled. 

In what follows we shall adopt the MV model as our initial condition at y — 0. The 
remaining part of this section is devoted to a study of the gluon distribution in this 
model, with emphasis on the Cronin peak, and, more generally, on the role of non-linear 
effects in rearranging the gluons in momentum space. Although this model has been 
extensively studied in the literature [10,12,14,29,30] (in particular, in relation with the 
Cronin peak [23-25, 53, 54]), our analysis below will bring some new results and conceptual 
clarifications, and will reveal some novel, and perhaps surprising, features, which to our 
knowledge have escaped to previous investigations. 

Specifically, in Sect. 2.2, we shall provide a complete, analytic, study of the gluon dis- 
tribution in the MV model, that we shall generalize on this occasion to include running 
coupling effects, in such a way to be compatible with the quantum evolution to be dis- 
cussed later. The results obtained in this analysis will permit us to clarify the conditions 
for the emergence of the peak, and explicitly compute properties like the location of the 
peak and its magnitude (in Sect. 2.3). Then, in Sect. 2.4, we shall reexamine a global 
argument in favor of the Cronin peak, based on a sum-rule [24] , that we shall use in order 
to better understand the origin of the gluons which make up the condensate. 

A rather surprising feature which will emerge from this analysis refers to the redistribution 
of gluons in momentum space under the influence of the non-linear effects : Whereas it 
has been since long appreciated that the effect of the repulsive interactions is to push 
the gluons towards the modes at larger momenta, and thus provide a spectrum which 
is infrared-safe [12,14,30], it has not been recognized so far that the dominant part of 
this spectrum — the one which is parametrically enhanced for large A, and provides a 
plateau of order l/a s at saturation (k± < Q S (A)) - - is actually compact, i.e., it falls 
off exponentially with increasing k± above Q S (A). In fact, as we shall see, it is precisely 
this exponential decrease in v?a(^_l) f° r momenta just above Q S {A) which explains why a 
pronounced peak appears in the ratio TZ p a when computed in the MV model. 

To conclude this general discussion, let us mention two alternative definitions for the 
unintegrated gluon distribution which are used in the literature. These definitions have 
in common the fact that they relate the gluon distribution to a scattering operator: the 
scattering amplitude Na(j'±,u) for a color dipole with transverse size r± which scatters 
of a hadronic target with atomic number A at relative rapidity y. Under suitable approx- 
imations, this quantity obeys a closed, non-linear, evolution equation (a generalization of 
the BFKL equation), originally derived by Kovchegov [33]. 

Specifically, in Ref. [38], the following definition has been proposed (the color dipole is 
taken to be made off two gluons): 
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a s N c 



(2.2) 



This definition has no deep motivation, but is simply based on the analogy with a formula 
for the true occupation factor (2.1) which holds within the MV model (see Sect. 2.2 
below). In fact, the quantities in Eqs. (2.2) and (2.1) are very similar to each other: They 
coincide with each other in the linear regime at high momenta k± ^> Q S (A, y), where they 
both obey the BFKL equation, and they also show a similar behavior at low momenta 
k± <C Q s (A,y), so they differ, at most, in the transition region towards saturation. The 
approximations that we shall develop in this paper are not sensitive to the details of 
this transition region, nor to the differences [47, 48] between the Kovchegov equation and 
the general, functional, evolution equation for the CGC, so all the results that we shall 
obtain apply literally to both definitions. Thus, our results are directly comparable to 
the numerical calculations in Ref. [25], which are based on Eq. (2.2) together with the 
Kovchegov equation. 

A different definition, with a deeper physical motivation, has been introduced in Ref. [19], 
and reads 



where in this context ipA(k±,y) is the function defined in Eq. (2.2), and SA(f±,y) = 
1 — ACi(r_L,y) is the ^-matrix element for dipole-hadron scattering. The quantity (2.3) 
enters a factorized formula for the cross-section for gluon production in proton-nucleus 
collisions [14,20-22], and as such it is directly relevant for the phenomenology of d-Au 
collisions at RHIC. Note that its interpretation as a 'gluon distribution' is only conven- 
tional (this is based on an analogy with the -factorization which holds in the linear 
regime at not so high energies), and should not give rise to confusion: There is a priori no 
reason why the (canonical) gluon distribution should enter the calculation of observables 
for high energy scattering. Indeed, in the high-density environment at high energy, the 
scattering operators are non-linear in the color field in the target (the non-linear effects 
describe multiple scattering), so, unlike what happens at low energy, they are not simply 
proportional to the 2-point function which defines the gluon distribution. 

The results that we shall obtain for ip A later in this paper do not directly apply to Ha, 
although they could be easily translated for it (at least, in the case of a fixed coupling), by 
using Ha oc k 2 V\(pA- The numerical analysis in Ref. [25], which considers both definitions 
(2.2) and (2.3), shows that the ratio 7Z p a is qualitatively the same when computed with 
any of these two definitions. But important quantitative differences persist, especially at 
low momenta, where the functions ipA and Ha are very different. In view of this, it would 
be interesting to repeat for h A) and also for the convolution yielding the cross-section for 
gluon production, the analysis that we shall give below in this paper for ipA- A potential 



h A (k±,y) = -k 2 ± Vl<p A (k±,y) = k\ 



' r -L c -ik ± -r ± 
47T 



a s N c 



(2.3) 
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difficulty that we foresee with such an analysis is the fact that the generalization of 
Eq. (2.3) to the case of a running coupling is ambiguous (unlike for the other definitions, 
Eqs. (2.1) and (2.2), for which a natural generalization exists; see below). 

2.2 The McLerran-Venugopalan model: Fixed & running coupling 

The color sources which generate the small-rr gluons in the McLerran-Venugopalan (MV) 
model [10,12] are the 3A valence quarks (from the A nucleons), which are assumed to 
be uncorrelated with each other except for the long-range correlations associated with 
confinement. The gluon occupation factor </?a(^_l) = <^A(k±,y = 0) at momenta k± ^> 
Aqcd is then obtained as [14, 30] 



where A is a non-perturbative scale of order Aqcd (the only trace of confinement), 
N c = 3 is the number of colors, and Q\ = a s N c fiA oc A 1 ^ is proportional to the color 
charge squared /j,a of the valence quarks per unit transverse area. The factor of 4 within 
ln(4/r^A 2 ) is only a matter of choice. The integration in Eq. (2.4) must be restricted to 
rj_ < 2/A, for consistency with the approximations leading to this formula [12], and also 
to avoid that the logarithm in the exponent changes sign. But as long as k± ^> A, the 
value of the integral is very little sensitive to the precise value of the upper cutoff. 

Since the MV model is a classical approximation, Eq. (2.4) is a priori written for a fixed 
coupling a s . Later, we shall consider quantum evolution with a running coupling, and at 
that stage we shall also need a generalization of Eq. (2.4) which includes running coupling 
effects to one-loop order (i.e., for a running coupling a s (Q 2 ) = b / \n(Q 2 / Aq CD )) . The 
generalization that we shall consider reads (from now on, we shall not distinguish between 
the scales A and Aqcd) 



where Q\ = boN c fiA has now a slightly different meaning as compared to the fixed coupling 
case, but is still proportional to A 1 / 3 . Formally, Eq. (2.5) is obtained from Eq. (2.4) by 
replacing in the latter a s — > a s (4/r[) within the denominator of the integrand, and also 
within Q\. In contrast to Eq. (2.4), there is no need for an upper cutoff in Eq. (2.5): the 
integral is convergent as written, and for any k± ^> A it is dominated by perturbative 
sizes rj_ <C 1/A. 




(2.4) 




(2.5) 
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Since the choice of a running in the absence of a complete one-loop quantum calculation 
is a priori ambiguous (especially in the presence of several scales, like in Eq. (2.4)), it is 
important to justify in more detail our choice leading to Eq. (2.5). Note first that (a) the 
coupling a s which is explicit in Eq. (2.4) (in the denominator) and the one which is implicit 
in the definition of Q\ have the same origin [12, 30], and thus need to be renormalized in 

the same way, and (b) the quantity 1— exp j — | r\Q\ In in the numerator of Eq. (2.4) 

can be recognized as the scattering amplitude for a color dipole of size r± which scatters of 
the nucleus. More precisely, this is the dipole made of the two primary gluons radiated by 
the valence quarks (the gluons which define the distribution), and the scattering amplitude 
describes the color precession of these primary gluons when propagating through the 
color field of the nucleus (see [12] for details). As well known, a small dipole couples 
predominantly to gluons with momenta k\ ^ 1/rL so it is indeed natural to choose the 
size of the dipole as the argument for the running of the corresponding coupling constant; 
this is also the choice made in other studies of the evolution of the dipole scattering 
amplitude with increasing energy [39,40,46,50,51]. 

There is one more subtle point about our choice in Eq. (2.5) : as a density of color charge, 
the quantity \ia involves itself a factor of a s , that we have implicitly treated as a constant 
in our arguments above. Specifically, for A x N c valence quarks which are homogeneously 
distributed within the nuclear disk with radius Ra, one obtains [12] 



= -52-- (2-6) 
K A 

Treating a s in Eq. (2.6) as a constant amounts to neglecting quantum corrections to the 
distribution of the valence quarks, which is in the spirit of the MV model. One sees that, 
in writing Eq. (2.5), we have treated differently the vertices describing the radiation of 
classical color fields from the valence quarks, and those describing the scattering of the 
color dipole off these classical fields. This is justified for the present purposes because these 
vertices play different roles in the subsequent evolution with increasing y: Whereas the 
valence quarks act simply as sources for the primary gluons at y — 0, and as such they 
can be viewed as classical colored particles, the quantum evolution proceeds via gluon 
radiation from the primary gluons, and the vertices describing this radiation are of the 
same type as those describing the scattering of the color dipole in Eq. (2.4). Thus dressing 
just these vertices is the minimal way to render Eq. (2.4) consistent with the quantum 
evolution which includes running coupling effects. 

The emergence of the dipole scattering amplitude in Eqs. (2.4) and (2.5) illustrates the 
deep connection between unitarization effects in scattering processes at high energy and 
saturation effects in the nuclear gluon distribution: Both types of effects arise from non- 
linearities associated with strong gluon fields, which here are encoded in the exponential 
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terms 2 in the integrands of Eqs. (2.4) and (2.5). The transition from the linear regime to 
the non-linear one takes takes place when the exponent is of order one. More precisely, we 
shall define the saturation momentum Q S (A) such that the exponent is equal to one when 
r = 2/Q s (A). For running coupling, Eq. (2.5) immediately implies Q 2 {A) = Q\ ~ A 1 / 3 . 
For fixed coupling, the saturation scale turns out to be larger than the scale Qa introduced 
by the color sources 



Q 2 S {A) = Q\ In 



A 2 



A l/3 In A 1 / 3 



(2.7) 



and this difference will be seen to have important consequences. Loosely speaking, the 
ratio 



Pa = 



In 



A 2 



In A 1 ' 3 



(2.8) 



will play the same role in the fixed coupling case as the inverse of the coupling constant 
in the case of a running coupling. In fact, for a running coupling, we shall define similarly 
Pa = hi(Q/i/A 2 ), an d then it is clear that 1/pa is essentially the same as the coupling 
evaluated at the saturation momentum : a s (Q 2 A ) = b /pA- Since the MV model makes 
sense only for weak coupling, it is justified to treat pa as a large parameter, pa ^> 1, 
which we shall do in what follows. It is interesting to note in this context that, for a 
gold nucleus at typical RHIC energies one expects Q 2 S (A) ~ 2 GeV 2 , which together with 
A ~ 200 MeV, yields p A ^ In 50 ~ 4. 



It is straightforward to compute the dominant behavior of the gluon distribution at asymp- 
totically large (k± ^> Q S (A)) or small (k± Q S (A), with k± ^> A though) momenta. 
One finds : 

(a) Fixed coupling: 



1 for k ± ^>Q s (A) 



a„N c k 



(a) Running coupling: 



(2.9) 



Physically, these exponentials represent the 5-matrix for dipole-nucleus scattering. 
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1 q2a for k ± ^>Q A 



The average inverse coupling constant which enters the last equation is defined as 



1 \ 1 (, Q\ M] k 2 ± \ 1 k ± Q A 



The high-momentum behavior in these equations, which is obtained after linearizing the 
exponential terms Eqs. (2.4) and (2.5), is recognized as the bremsstrahlung spectrum due 
to radiation from independent color sources. This would be the spectrum at any k± in 
the absence of non-linear effects in the dynamics of the radiated gluons. But already at 
high-fcx, this spectrum receives non-linear corrections which are suppressed by powers 
of Q 2 s {A)/k\ ('higher-twist' effects; see below). The dominant behavior at low momenta 
is obtained after neglecting the exponential terms Eqs. (2.4) and (2.5), and is correct to 
leading logarithmic accuracy; that is, the first corrective term would be a constant term 
under the logarithm. 

As we shall see in Sect. 2.3, the asymptotic expressions above are already sufficient to 
demonstrate the existence of a peak in the ratio 7Z p a- But in order to study the properties 
of this peak, we also need the gluon occupation factor at generic (intermediate) momenta, 
which given the simplicity of the MV model can be computed analytically. 

In fact, in the case of a running coupling, we have been able to evaluate the integral in 
Eq. (2.5) in analytic form, with a result which is displayed in the Appendix, and which 
will be also discussed below. For fixed coupling, we do not have such an exact result, 
but we have found a convenient series representation of the integral in Eq. (2.4), which 
converges very fast, and renders the physical interpretation of the result transparent. 

Let us consider the fixed coupling case first. After rewriting the exponent in the integrand 
of Eq. (2.4) as (cf. Eq. (2.7) and (2.8)) 



one can decompose the integral into two pieces as 
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Pa , 

= <ftf(k ± ) + tpT k ik±), (2.13) 
where the first, 'saturating', piece can be explicitly evaluated, and reads 



<ftf(k ± ) = \t r(0, z) , (2.14) 

a s N c 

while the second piece - - the sum of all-'twist' contributions (see below) - - can be 
evaluated as a series expansion, obtained by expanding the exponential within the square 
brackets in powers of (tlnl/t)/p j4 as 



1 r , ™ max / t \ 11 1 In™ t 

<pT\k ± ) = — / dtM^e-* E — —■ (2-15) 

a s N c p A J n^i\PAj n\ 

In these equations, z = k\/Q 2 8 (A), T(0,z) is the incomplete Gamma function: 



POO p t 

r(0,z)= —dt, (2.16) 

J z t 

n mffl ~ e PA = Q 2 s (A)l 'A 2 ^> 1, and the truncation of the series in Eq. (2.15) at n < n max 
reflects the upper cutoff r max = 2/ A which is implicit in the integral over rj_ in the 
first line of Eq. (2.13). (The corresponding integral giving iffi, Eq. (2.14), has been 
extended up to infinity, since rapidly convergent.) Without a truncation, the series would 
be asymptotically divergent, but the divergent behavior would start to manifest itself only 
at very large orders, beyond n max . On the other hand, when p A 3> 1 the truncated series 
is rapidly convergent for any z, and therefore ^(^jj can be computed even analytically 
with high accuracy (see the Appendix for details). 

The first important observation about the decomposition in Eq. (2.13) is that the second 
piece, v^a" 8 *) is multiplied by an overall factor 1/pa (this is manifest on Eq. (2.15)), and 
thus is parametrically suppressed at large A. This implies that, in the saturation region 
at z ^ 1, the gluon distribution is dominated by the first piece, ip s A . In particular, 
(p^ captures the leading behavior at low momenta k± <^ Q S (A), cf. Eq. (2.9): indeed, 
r(0, z) ~ \nl/z for z< 1, whereas cp^ 1 ^ is analytic near z = 0. In fact, if one remembers 
that 1/pa can be effectively identified with a s (cf. the discussion after Eq. (2.8)), it 
becomes clear that the l/o; s -enhancement of the gluon distribution, which is the hallmark 
of saturation, is associated solely with (p^ : For z < 1, the latter provides a plateau with 
a height of 0(l/a s ), whereas y^™ 3 * brings only a small correction, of 0(1). 
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The second important observation is that the saturating piece yf^ is compact, in the sense 
that it vanishes exponentially at momenta outside the saturation region (since T(0,z) « 
e~ z /z for z >> 1). By contrast, it can be checked that, for z 1, the function ip^ 1 ^ can 
be expanded in powers oil/ z (up to logarithmic corrections), thus generating the 'twist 
expansion' of the MV model. It should be stressed that each of the terms in the truncated 
series in Eq. (2.15) contains terms to all orders in the twist expansion: the first term 
(n = 1) includes the bremsstrahlung spectrum oc 1/z, together with an infinite series of 
'higher-twist' terms (i.e., terms of order 1/z 2 , or higher), the second term (n = 2) starts at 
order 1/z 2 , etc. The first few terms in the twist expansion are exhibited in the Appendix. 

We see that the gluon distribution of a large nucleus in the MV model naturally decouples 
into two pieces, one which controls the spectrum at saturation, and another one which 
controls the tail of the distribution at high momenta. It is natural to interpret the first 
piece, yr^, as the occupation factor in the color glass condensate. The fact that this dis- 
tribution appears to be compact is probably just specific to the MV model, since, as we 
shall see, this feature is washed out by the quantum evolution (cf. Sect. 5.3). The previous 
considerations also suggest that the twist terms contribute only little to the physics of 
the Cronin peak in the MV model. This will be confirmed by the analysis in the next 
subsection. 

The corresponding discussion in the case of a running coupling is even simpler, since the 
corresponding 'twist' piece involves only one term, instead of a series of terms. Also, some 
arguments become more transparent since, with a running coupling, l/a s (Q 2 (A)) and 
Pa/% are explicitly identified. Specifically, Eq. (2.5) can be decomposed as: 



where the saturating piece is explicitly enhanced by the factor p A (that is, it is of 0{l/a s )), 
while (/9^ lst is of 0(1) for z ~ 1, and can be expanded in powers of 1/z for z 3> 1 (thus 
generating the twist expansion). The only noticeable difference with respect to the fixed 
coupling case is that, now, ip t ^ mt {z) is not analytic near z — > 0, but rather yields a large 
negative contribution —{l/2bQN c )\Yi'z in that limit, which contributes to the dominant 
low-momentum behavior exhibited in Eq. (2.10). Since, however, the overall contribution 
remains positive for any k± > A, it is clear that (p s £ t (k±) is still the dominant contribution 
at any A < k± < Qa- The integral giving < ! p^ vlst is explicitly evaluated in the Appendix, 
where its twist expansion will be also considered. 

The previous considerations are illustrated in Figs. 2 and 3, which exhibit the gluon 
occupation factor </?a(&±), and its various contributions y^f^/cjj and (p^ st (k±), m the 
MV model with fixed and running coupling, in either linear scale (Fig. 2), or log-log scale 




(2.17) 
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in a wider range of momenta (Fig. 3). 




z 



12 3 4 

9A 




z 



12 3 4 



Fig. 2. The gluon occupation factor <pa{z) as a function of the scaled momentum variable 
z = k 2 /Q 2 s (A) in the MV model with either hxed (figure above) or running (below) coupling 
and pa = 6. With respect to the text definitions, we plot the quantities pAOt s N c x ifA for fixed 
coupling, and boN c x ipA for running coupling. The black (thick) line corresponds to (pa(z); the 
blue (solid) line shows the saturation contribution (p^^z); the magenta (dotted) line shows the 
twist contribution ip A Mlst (z); the green (dashed) line represents the bremsstrahlung spectrum. 
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Fig. 3. Logarithmic plot of the gluon occupation factor <pa{z) as a function of the scaled momen- 
tum variable z = k 2 /Ql(A) in the MV model with either fixed (figure above) or running (below) 
coupling and pa = 6. The plotted quantities are rescaled as explained in the caption to Fig. 2. 
The black (thick) line corresponds to the gluon occupation factor (pA_(z); the blue (solid) line 
shows the saturation contribution (p 3 ^ (z); the green (dashed) line represents the bremsstrahlung 
spectrum. 
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2.3 Cronin effect in the McLerran-Venugopalan model 



To form the ratio TZ p a, Eq. (1.1), one also needs the gluon distribution in the proton, 
ip p (k±,y). For y = and k± ^> A, the proton is in the perturbative (i.e., linear) region, 
and will be described by the bremsstrahlung spectrum in Eqs. (2.9) or (2.10) (for fixed 
and running coupling, respectively) with Q\ — > Q 2 , and Q p = 0(A). For instance, 



1 Q 2 

<P P (h±) = — — -pf 6(/cj_ — Q p ) (fixed coupling). (2.18) 

Oi s l\ c Kj_ 

We shall assume the following relation between the two scales: 



Q\ = A^Ql (2.19) 

which holds formally if one extrapolates the MV model (a priori valid for a large nucleus) 
down to A = 1. 



Note that the quantity in the denominator of Eq. (1.1), namely A 1 ^ 3 (p p (k±), coincides 
with the nuclear gluon distribution in the high momentum limit (i.e., the bremsstrahlung 
spectrum in Eqs. (2.9) or (2.10)). Thus, within the MV model at least, the ratio TZ p A{k±) 
is also a measure of the deviation of the actual nuclear gluon distribution from the cor- 
responding prediction of linear perturbation theory. As already discussed, this deviation 
is associated with non-linear effects in the gluon dynamics, in particular with saturation. 
We immediately conclude that lZ p A(k±) must approach one at high k±: 



K pA (k ± ) ~ 1 for k ± > Q S {A). (2.20) 

By using this condition together with very general properties of saturation 3 , we can 
deduce the existence of a Cronin peak without any detailed calculation: 



Consider the fixed coupling case, for definiteness. The basic consequence of saturation is 
that, for momenta k± ^ Q S (A), the nuclear occupation factor develops a plateau with 
the height of 0(l/a s ); that is, v?a(^_l) ~ {^-/ ot s N c )^>{z), where z = k\/Q 2 s {A), and the 
function (p(z) is slowly varying when z < 1, and is of 0(1) when z ~ 1. This function 
depends upon k± and A only through z because it is dimensionless, and Q S (A) is the only 

3 Of course, all the properties to be discussed below are manifest on the expressions (2.14) and 
(2.17) for (f S A t (k±), but here we would like to construct our argument without relying on the 
explicit formulae for ipA that we have obtained previously, in order to emphasize the general 
conditions required by the existence of the peak. 
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scale in the problem other than k± (since at saturation the spectrum cannot be sensitive 
to the non-perturbative scale A). We deduce that: 



K pA {k±) « 7^-^(3) = PAZ<p(z), for k ± <Q s (A). (2.21) 

Since (p(z) is only logarithmically divergent as z — > 0, it is clear that the ratio lZ p A(k±) is 
much smaller than one for k± small enough (e.g., it is of 0(A" 1 / 3 ) when k± ~ Q p ), but 
it is of 0{pa) — and thus strictly larger than one — for k± = Q S (A). This behavior at 
small k±, together with the asymptotic behavior (2.20) at large k±, immediately imply 
that the ratio 1Z p a must have a maximum at some intermediate value of k±. 

This general argument does not tell us where is the maximum actually located, and not 
even whether there is a single maximum, or several. But it is highly probable, and will 
be confirmed by the explicit calculations below, that there is indeed only one maximum, 
which is located near Q S (A) and thus has a height of O(pa) '■ Indeed, there is no intrinsic 
scale larger than Q S (A) in the problem, and for k± < Q S {A) the ratio is still increasing 
with k±, as manifest on Eq. (2.21). 

Note the specific way how the large factor pa has entered the calculation in Eq. (2.21) : 
this is due to the mismatch (2.7) between the scale Q\ associated with the color sources 
and the saturation scale Q 2 S {A) generated by the non-linear gluon dynamics. It is easy 
to check that a similar enhancement occurs also for a running coupling, although in that 
case the argument is more direct: the factor pa is introduced by the inverse coupling in 
the saturation condition (pA(k± = Qa) ~ ^/ a s (Q 2 A )N c , with k ~ 0(1). 

We conclude that the existence of the Cronin peak is a consequence of the fact that the 
gluon distribution at saturation is larger, by a factor of pa ~ In A 1 / 3 , than the naive 
extrapolation of the bremsstrahlung spectrum down to kj_ ~ Q S (A). This logarithmic 
enhancement is due to non-linear effects which cause the gluon occupation factor to 
increase faster with l/k\ than predicted by linear perturbation theory, before eventually 
saturating to a value of order l/a s at k± ^ Q S (A). 

We now rely on the explicit formulae for tfA established in the previous subsection to 
give a complete description of the peak. For more clarity, we start with the large- A 
limit, in which the gluon occupation factor for momenta around Q S (A) is given by the 
saturating piece (p^ 1 , up to corrections of 0(1/ pa)- Using (p A t (k±) from either Eq. (2.14) 
or Eq. (2.17), one finds the same expression for the ratio 1Z p a with both fixed and running 
coupling, namely: 

n pA (k ± ) w zT(0,z)p A for k± ~ Q S (A) . (2.22) 
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This is of the form anticipated in Eq. (2.21), and thus is of order pa when z ~ 1 (T(0, 1) = 
0.219...). It can be checked that the function f(z) = zF(0,z) has a maximum at z « 
0.435, with the peak value f{zo) ~ 0.281. Moreover, this peak is very pronounced, since 
f(z) is almost linearly increasing below the peak, but exponentially decreasing above it. 

It is then straightforward to add the contribution of the twist terms in Eq. (2.15) or (2.17), 
and compute their effect on the location and the magnitude of the peak in an expansion 
in powers of 1/ Pa- One finds (with 7^ max (^4) = 71 p a(zo)) : 



0.882 2.122 / ,x 

z = 0.435 + + — =- + 0[pa 3 ), 

Pa Pa v 1 

294 / \ 

7^max(^4) = 0.281 pa + 0.300 + + 0(p A 2 ), (2.23) 

Pa v 7 

in fixed coupling case, and respectively 



1.382 2.038 / ox 

z = 0.435 + + — =- + 0[pa 3 ), 

Pa Pa V ' 

TWA) = 0.281 p A + 0.524 + — + 0(p A 2 ), (2.24) 

Pa v ' 

in the running coupling case. 

Finally, for both fixed and running coupling one can compare the results obtained above 
in the 1/pa expansion to the corresponding exact results, and the agreement is good even 
for pa = 4 (which, we recall, is a realistic value at RHIC). For example, with running 
coupling and pa = 4, the exact calculation yields z = 0.89 and 7£ max = 1.85, while 
Eq. (2.24) predicts z = 0.91 and ft max = 1.85. 



In Fig. 4, where we show the Cronin ratio in the MV model with both fixed and running 
coupling, we also represent the individual contributions separated in the r.h.s. of Eq. (2.13) 
(or (2.17)), to better emphasize that the 'twist' effects give only a small contribution in 
the region of the peak. As obvious on this figure, the effect of the twist corrections is to 
flatten the peak a little bit, and also to move it towards larger momenta. Still, a well 
pronounced peak emerges because, even for pa as small as 4, the dominant contribution 
to v?a(^_l) f° r k± around Q S (A) is still given by the compact, 'saturating', piece (p S A- 
Accordingly, the nuclear occupation factor shows an exponential fall-off for momenta just 
above Q S (A). It can be easily checked that the behavior changes from an exponential to 
a power law around z ~ lnp^- 
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Fig. 4. The Cronin ratio 1Z p a{z) as a function of the scaled momentum variable z = k 2 /Q^(A) 
in the fixed (above) and running (below) coupling McLerran-Venugopalan model for pa = 6. 
The black (thick) line corresponds the ratio 1Z p a(z); the blue (solid) line shows the saturation 
contribution IZ^(z); the magenta (dotted) line shows the twist contribution TZ p w A mt {z). 
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2.4 A sum rule and its consequences 



In this subsection we shall discuss an alternative, global, argument in favor of the existence 
of the Cronin peak in the MV model, which sheds more light on the role of non-linear 
effects in the nuclear gluon distribution. This argument, due to Kharzeev, Kovchegov, and 
Tuchin [24], is specific to the MV model — it is based on a sum-rule which reflects the basic 
property of this model that color sources are uncorrelated — , and cannot be extended to 
y > (since the quantum evolution introduces correlations among the color sources; see 
below). But even within the MV model, the use of this argument is quite subtle, because of 
possible complications with ultraviolet divergences. In fact, as we shall explain below, the 
original formulation of the relevant sum-rule in Ref. [24] is not mathematically rigorous, 
which may have given rise to confusion (especially, in relation with the generalization of 
this sum-rule to y > 0; see the discussion in Ref. [25]). In what follows, we shall properly 
restate the argument, and then present a stronger version of it, based on partial sum-rules, 
which provides a more detailed information on the distribution of gluon at saturation. 

Specifically, the argument in Ref. [24] is based on the following 'sum-rule' (the IR cutoff 
A is strictly needed only for the proton) : 



which holds indeed within the MV model, for both fixed and running coupling, as we 
shall demonstrate shortly. From this relation, one can infer the existence of a Cronin peak 
through the following reasoning: Since at low momenta k± <C Q S (A) we have <^a(^±) < 
A 1 / 3 (p p (k±) due to gluon saturation in the nucleus, whereas at large momenta k± Q S (A) 
the two functions approach the same limit, one concludes that <fA(k±) must be larger than 
A 1 / 3 ip p (h±) at some intermediate momenta, in order to give rise to the same integrated 
distribution. One reason for the debate around this argument in the literature [25] is that 
in Ref. [24] the sum-rule has been abusively written as: 



At a first sight, this might look identical to Eq. (2.25), but in reality Eq. (2.26) is mean- 
ingless, since the expressions on both sides are logarithmically divergent in the ultraviolet. 
If, moreover, one introduces an upper cutoff Q 2 to eliminate these divergences, then the 
ensuing, finite, expressions do not coincide with each other for any finite Q 2 (see below). 
The would-be divergent pieces (i.e., the contributions proportional to InQ 2 ) are indeed 
the same for both integrals — because (fA(k±) — A 1 / 3 ip p (k±) for large k± Q 3 (A) — , 
but by itself this property only tells us that the integral in the left hand side of Eq. (2.25) 




(2.25) 




(2.26) 
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is well denned, but not also that the value of this integral is actually zero. In other terms, 
the equality of the divergent pieces in Eq. (2.26) carries no information about the behavior 
of the integrands at intermediate momenta, and thus cannot guarantee the existence of 
the Cronin peak [25] . 

Still, as we shall prove below, the difference Ga(Q 2 ) — A l l z Q p (Q 2 ) vanishes like l/Q 2 when 
Q 2 — > oo, so the sum-rule (2.25) holds indeed as written. We have defined here: 



rQ 2 d 2 k , 

g A (Q 2 ) = / -^L<p A (k ± ), (2.27) 

J A 2 71 

(together with a similar expression for the proton) which, up to a global factor tcR a x 
(N 2 — l)/(27r) 2 (cf. Eq. (2.1)), is the integrated gluon distribution in the MV model, i.e., 
the total number of gluons localized within a transverse area l/Q 2 - 

Our proof of Eq. (2.25) is very similar to the original one in Ref. [24], namely it uses 
the Fourier representation of the integrand (cf. Eqs. (2.4)-(2.5)) to easily perform the 
integral over k±, and thus produce a 5-function 5^ 2 )(rj_). Then, the sum-rule immediately 
follows because the integrand vanishes when r± — > 0. For instance, in the case of a running 
coupling, Eq. (2.5), we write (we omit unnecessary constant factors) : 




(2.28) 



which vanishes as r\ lnr^ when r± — > 0. 

If the above integral over k± is restricted to some finite, but large Q 2 , with Q 2 3> Q\, 
then instead of the 5-function 5 ( - 2 \r±) one generates a smooth distribution in r± which is 
peaked at rj_ = 0, with the peak height proportional to Q 2 , and the width of the peak of or- 
der l/Q 2 . When integrated over such a distribution, the 'higher-twist' terms in Eq. (2.28) 
(i.e., the higher order terms in the expansion of the exponential) give contributions which 
vanish as powers of l/Q 2 when Q 2 — > oo. The dominant higher-twist contribution, of 
order l/Q 2 , turns out to be negative (this is explicitly computed in the Appendix). Thus, 
for a finite, but large, Q 2 , the total number of gluons with transverse area l/Q 2 is smaller 
in the MV spectrum of a nucleus than in the corresponding bremsstrahlung spectrum 
(i.e., than it would be in the absence of non-linear effects). Specifically, one can use the 
results in the Appendix to deduce that 
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Qa{Q 2 ) - Qbs{Q 2 ) 



1 01(A) 



for Q 2 ^Q 2 (A), 



(2.29) 



b N c 2Q 2 



with Qbs{Q 2 ) = A l ' 3 Q p {Q 2 ). This is illustrated in Fig. 5. 

The sum-rule (2.25) can be understood as the consequence of two basic properties, out of 
which one is generic — the fact that the non-linear effects fall off as inverse powers of k\ in 
the high-momentum limit — , and the other one is specific to the MV model: the fact that 
the color sources (the valence quarks) are uncorrelated. Specifically, the second property 
implies that, at high momenta k± ^> Q S (A), the gluon distributions produced by the 3 A 
valence quarks sum up incoherently, so that the leading-twist terms cancel exactly in 
the difference (fA(k±) — A 1 / 3 ip p (k±). Moreover, by the first property, the remaining terms 
in this difference involve higher powers of which explains Eq. (2.29), and thus the 

sum-rule (2.25). 

The above considerations also help us to understand why the sum rule is bound to fail 
after taking the quantum evolution with y into account: At y > 0, the color sources are 
predominantly gluons which are themselves products of radiation, and thus are correlated 
with each other (unlike the valence quarks) even in the leading-twist approximation. 
Because of that, the gluon occupation factors in the nucleus and in the proton are not 
simply proportional to each other, not even at high k±. This can be easily verified on 
Eqs. (3.15) and (3.26), which give the dominant behavior at large momenta for y > 0. For 
instance, for fixed coupling, Eq. (3.15) implies (we omit the trivial factor l/a s N c ) : 



where p(A, k±) = In k\/Q 2 s {A), p{p, k±) = \nkj_/Q 2 = p(A, k±) + pa, and in writing the 
second line we have kept only the leading term at large p(A, k±). Although the expression 
above does vanish as k± — > oo, it is clear that it is not integrable in the sense of Eq. (2.25). 
The very fact that the sum-rule fails to apply at y > does not preclude the existence 
of the Cronin peak: The sum-rule is only a sufficient condition for the existence of the 
peak, but not also a necessary one. And indeed the peak does persist (although its height 
is rapidly decreasing) in the early stages of the evolution, as we shall show in Sect. 5. This 
is also seen in the numerical results in Ref. [25]. 

Returning to the MV model, where the sum-rule (2.25) is valid, it is instructive to un- 
derstand its physical implications in more detail: Eq. (2.25) tells us that the net result 
of the non-linear effects in the gluon dynamics is merely a redistribution of the gluons in 
the transverse momentum space. As it should be clear from the previous arguments, and 



(fiA(k±,y) - A 1/3 (p p (k±,y) 




(2.30) 
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Fig. 5. The integrated gluon distribution function Qa(Z) as a function of the scaled mo- 
mentum variable Z = Q 2 /Q 2 S (A) in fixed (above) and running (below) coupling McLer- 
ran-Venugopalan model for pa = 6. With respect to the text definitions, we plot the quantities 
\pAOt s N c / Ql(A)] x Qa for fixed coupling, and [boN c /Q 2 (A)] x Qa for running coupling. The black 
(thick) line corresponds to the integrated gluon distribution function Qa(Z); the blue (solid) line 
shows the saturation contribution Q s ^ t (Z); the green (dashed) line represents the bremsstrahlung 
distribution function. 
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also from the plot in Fig. 3, this redistribution is associated with repulsive interactions 
which push the gluons towards the modes with larger momenta. In Fig. 3 one can see 
that, whereas at low momenta k± < Q C (A) the gluon occupation factor is larger for the 
bremsstrahlung spectrum than for the MV spectrum, at higher momenta k± > Q C (A) the 
opposite situation occurs. Here, Q C (A) is defined by 4 (see Eq. (2.9)): 



^~ln%^ for k^Q c {A\ (2-31) 
which implies : 



2 (A) ~ ® A Qa „ Q 2 s( A ) (2^2) 



'A 



where we have also used Q 2 S {A) = Q a Pa, cf. Eq. (2.7). It is easily checked that A 2 <C 
Q 2 C (A) < Q 2 S (A) (since p A > 1). 

Thus, the effect of the interactions is to move the gluons up in transverse momentum, 
from below Q C (A) to above it. This explains, in particular, why when counting the gluons 
up to a finite (large) Q 2 , one necessarily finds a smaller number for the MV spectrum 
than for the bremsstrahlung one, cf. Eq. (2.29). But as we shall explain now, at large 
k± ^> Qs(A) the effects of this redistribution are only tiny: Most of the gluons which are 
in excess in the bremsstrahlung spectrum at k± < Q C (A) are found in the MV spectrum 
in the saturation region at k± ^ Q S (A). In particular, the enhancement factor l/a s ~ pa 
characteristic of saturation finds its origin in the phase-space for bremsstrahlung radiation 
at momenta A < k± < Q C (A). This factor is large because of the infrared sensitivity of 
the bremsstrahlung spectrum. 

A simple way to see this is to notice that the total number of gluons in the region 
A < k± < Q S (A) is essentially the same for the bremsstrahlung spectrum and for the MV 
spectrum. Indeed, we have first : 



Qbs{Q s {A)) - J a2 —^ c1 ^-^- c ln ^r- - > (2-33) 

where the saturation momentum Q S (A) has been reconstructed as Q 2 S {A) = Q\pa, cf. 
Eq. (2.7). Note that the large factor pa = \nQ 2 (A)/A 2 has been generated from the 



4 The subsequent formulae in this subsection are written for the case of a fixed coupling, for 
definiteness. It can be easily checked that analog formulae, leading to similar conclusions, apply 
also for the MV model with a running coupling. 
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phase-space for the bremsstrahlung gluons at A < k± < Q S (A). For the MV spectrum, 
we shall define Qa{Q 2 s {A)) as the total number of gluons contained in the saturating 
component of the spectrum (p s ^ t (k±). This is appropriate since ifiA(k±) ~ t P B A t (k±) f° r an y 
k± < Q S (A), whereas at larger momenta <p A ^(k±) falls off exponentially. One then obtains: 



g Am = / ^ -L r(0, ,1/QKA)) = 2j£ . (2.34) 

(The simplest way to obtain this result is to notice that the unrestricted integral of 
tp s X t {k±) is tantamount to letting rj_ — > in the integral representation of (/^(/cjj i n 
Eq. (2.13).) As anticipated, the final expressions in Eqs. (2.33) and (2.34) are identical. 
But, of course, the way how these gluons are distributed in k± is very different in the two 
cases: For the bremsstrahlung spectrum, which is infrared divergent, most of these gluons 
are located at low momenta k± < Q C (A), while for the MV spectrum they are rather 
uniformly distributed at all momenta up to Q S (A). For instance (see Eq. (2.32)): 



GssiQKA)) = 1 _ A = _^. ]n _sLi w _^. {pA _ lnpA) , (2 .35) 



is almost the same as the result in Eq. (2.33) (in particular, it includes already the 
enhancement factor p& due to the logarithmic phase-space 5 ). This confirms the fact that 
the color glass condensate in the MV model is formed by redistributing the 'infrared' 
gluons at k± < Q C (A) in the bremsstrahlung spectrum. 



3 Non linear gluon evolution in the Color Glass Condensate 



Within the effective theory for the CGC [26], the evolution of the gluon distribution 
with y is not described by a single, closed, equation, but rather by an infinite hierar- 
chy of coupled equations, which relate n-point functions with arbitrary n > 2 (by itself, 
(p{k±,y) corresponds to a 2-point function), and are compactly summarized in a func- 
tional evolution equation, the JIMWLK equation. Fortunately, this complicated evolution 
simplifies considerably in interesting limiting cases, which are tractable through analytic 
approximations [35,36]. 

At low momenta k± ^ Q s (A,y), the gluons form a condensate with occupation factors 
of order l/a s , and the dynamics is highly non-linear. Still, one can rely on mean field 

5 Although Ql(A) « Q 2 s {A)/pa < Ql(A), the logarithmic phase-space at A < k ± < Q C {A) 
(namely pa — ^Pa) is almost the same as that for the whole interval A < k± < Q S (A) (which 
is equal to pa)- 
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approximations to deduce a formula for tp(k±,y) valid in the low momentum limit k± <C 
Q s (A,y) [35-37,55,58], which will be presented in Sect. 3.1 below. 

At high momenta k± ^> Q S (A, y), the gluon density is low, and non-linear effects in the 
evolution become negligible. Then, the general equations for the evolution of the gluon 
distribution boil down to a single, linear, equation: the BFKL equation [43]. The latter 
controls also the approach towards saturation from the above (at least, approximately), 
and thus determines the saturation momentum Q S (A, y) [7,38-41]. 

For a fixed coupling, we shall discuss in Sect. 3.2 the solution to the BFKL equation 
in the saddle point approximation and with saturation boundary conditions [38-41] (see 
also [48]). Our derivation of this solution will be quite schematic, and only intended 
to emphasize a subtle point — namely, the fact that it is the saturation scale in the 
initial conditions, Q S (A), which sets the infrared cutoff for the transverse phase-space 
available for evolution [42] — , and to prepare the ground for a general discussion of 
high-A;^ suppression through quantum evolution, to be presented in Sect. 4. 

For a running coupling, with the scale for running set by the gluon transverse momentum 
k±, we shall follow the analysis in Refs. [39-41], from which we shall simply quote here 
the relevant results, and discuss their range of validity. 

3.1 Non-linear evolution at low k± : Saturation and universality 

In the high density regime deeply at saturation (k± <C Q s (A,y)), the color fields are 
strong (A 1 ~ 1/g) and the dynamics is fully non-linear, which makes it difficult to use 
standard perturbative techniques. On the other hand, these conditions are favorable to 
the use of mean field approximations [35,36], which enable us to determine the y— and 
k±— dependencies of the gluon occupation factor. Similar results are also obtained from 
studies of the unitarization effects in dipole-hadron scattering [37,55,58]. 

Remarkably, it turns out that the functional form of the gluon occupation factor at sat- 
uration is universal [35]: For k± <C Q s (A,y), (pA{k±,y) is independent of either the initial 
condition at y — 0, or the details of the evolution leading to saturation, except for the 
corresponding dependencies of the saturation scale itself. Moreover, this universal form 
happens to be the same as in the MV model (cf. the expressions in the second lines of 
Eqs. (2.9) and (2.10)), which is quite non-trivial, and to some extent surprising, since 
the physical conditions leading to this functional form are very different in the two cases: 
Whereas in the MV model the color sources are uncorrelated, and the logarithmic behav- 
ior at low kj_ is merely the result of non-linear effects in the classical equations of motions 
for the color fields [12, 30], in the effective theory at small x the color sources are strongly 
correlated, in such a way to ensure color neutrality over a transverse size ~ l/Q s (A,y) 
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[35,36,58], and the logarithmic behavior at low k± is already encoded in the correlations 
of the color sources (i.e., it holds independently of the presence of non-linear terms in 
the classical field equations; the latter affect only the overall normalization of the gluon 
distribution, which in this non-linear regime is anyway not under control). 

Specifically, for fixed coupling, one finds [35-37] : 



<PA(k±,y) a* -^7- hi ^ 8 ^'^ , (fixed coupling), (3.1) 

a s M c kj_ 

whereas for running coupling (with the scale for running set by the gluon momentum k±) 
one obtains [12, 58]. 



Mk±,y)^^T (^) ln ^^' ( running cou P lin s)> (3-2) 

with (compare to Eq. (2.11)) : 



W^l^+^l' (3 - 3) 

The overall factor ao in these equations is a number of order one which cannot be computed 
analytically (and which needs not be the same for fixed and running coupling). 

Note that, in the fixed coupling case, the gluon distribution at saturation depends upon 
the two kinematical variables kj_ and y only via the ratio z = k\/Q 2 s {A, y). This property, 
known as geometric scaling [56], reflects the fact that the saturation momentum is the 
only intrinsic scale at saturation. As manifest on Eq. (3.2), the running of the coupling 
introduces a second intrinsic scale, namely A, and this is a source of geometric scaling 
violations, which are however under control. 



3.2 Linear evolution at high kj_ : Fixed coupling 



In the low density regime at high momenta k± ^> Q s (A,y), the dynamics is linear, and 
the nuclear gluon occupation factor can be obtained (to the present accuracy) by solving 
the BFKL equation [43] with initial conditions at y = provided by the MV model (cf. 
Eq. (2.4)-(2.5)) and with an 'absorptive' boundary condition at k± ~ Q s (A,y) [39,40] 
which mimics the non-linear effects in the approach towards saturation. (Similar results 
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have been recently obtained [41] through direct studies of the non-linear Kovchegov equa- 
tion [32,33].) In particular, the saturation scale itself is obtained by requiring (fiA{k±,y) 
to become of order l/a s at saturation (up to subleading terms of 0(1)) : 



<PA(k±,y) 



a s (Ql(A )y ))N c 



for k ± rsj Q s (A,y) . 



(3.4) 



(k is a number of order one.) 

When y and/or k± are relatively large (see below), an approximate analytic solution can 
be obtained, via a saddle point approximation. Here we shall outline only a few steps in 
the construction of this solution (see Refs. [38,39] for details). 

For fixed coupling, the solution to the BFKL equation for <fA(k±,y) can be written in 
Mellin form as: 



where a s = a s N c /n, x{l) = 2^(1) — VKt) — ^(1 ~~ 7) "with VKt) = d\nF('j)/d , y and 
complex 7 is the eigenvalue of the BFKL kernel, <^a(7; Qo) is the Mellin transform of the 
initial condition <£a(&±), Eq. (2.4): 



and Qo is an arbitrary reference scale introduced for dimensional reasons, and which drops 
out in the complete result, as obvious on the above equations. Given the behavior of the 
initial distribution (pA(k±) in various limits, cf. Eq. (2.9), it can be checked that the integral 
over kj_ in Eq. (3.6) is absolutely convergent for < Re 7 < 1. Thus the contour for the 
complex integration in Eq. (3.5) can be chosen as C = {7 = 71 + 272 ; —00 < 72 < 00}, 
with a fixed 71 in the range < 71 < 1. 

The same study of the convergence properties of the integral in Eq. (3.6) tells us that this 
integral is dominated by momenta of the order of the initial saturation momentum Q S (A): 
Indeed, for large k± ^> Q S (A), <fA{k±) oc l/k\, and the integral is saturated by momenta 
of the order of the lower cutoff Q S (A). Furthermore, for low momenta k± Q S (A), the 
initial distribution saturates: (pA(k±) oc \n(Q^/k\), so the integral is now dominated by 
momenta near the upper cutoff, i.e., by k± ~ Qs(A) once again. Thus, after perform- 
ing the integral over k±, Q S (A) replaces Qo as the natural reference scale in the Mellin 
representation of the solution ifA(k±,y). Then, Eq. (3.5) can be rewritten as: 




e a -«*M<^( 7 ;Q ), 



(3.5) 




(3.6) 




<Pa(t,Pa), 



(3.7) 
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where any reference to the arbitrary scale Qq has disappeared. Note that, for more 
clarity, we have extracted a factor l/a s N c out of the (Mellin transform of) the initial 
condition ^4(7; Pa)- The latter still depends upon Q S (A), but only logarithmically, via 
p A = \n(Q 2 s (A)/ A 2 ). The precise form of the function ^4(7; Pa) will not be needed in 
what follows. Rather, we shall use its approximate form which is obtained when </?a(&±) 
in Eq. (3.6) is replaced by the simplest interpolation between the asymptotic behaviors 
shown in Eq. (2.9). This reads : 



r(7 Mi-7) + £(7) 

Pa 7 Pa 1 - 7 ' 7 2 



^(7; pa) = - tiniV-JU + till + . . . „ + J_ . (3.8) 



The first term in the last (approximate) equality, which has a pole at 7 = 1 and is 
suppressed as 1/pa, comes from the high-fcj_ behavior (the bremsstrahlung spectrum in 
Eqs. (2.9)), while the second term, with a double pole at 7 = 0, comes from the saturating 
behavior at low k 1 . 



3.2.1 The saddle point approximation 

At this stage, it becomes natural to introduce the logarithmic momentum variable p = 
\n{k\/Q 2 s {A)) = p(A,k±), and notice that, when either d s y, or p, or both, are large, the 
integral in Eq. (3.7) can be evaluated via a saddle point approximation. Specifically, if one 
writes: 



VM-Zkl&^m™), (3-9) 



where: 



F(l,r) = -7r + x(7), 
then the saddle point 70 = 7o(r) satisfies the condition: 

= -r + x'(7o(r)) = 0, 



(3.10) 



dF 



(3.11) 



7=70 (r) 



which shows that 70 (r) is a real number, in between and 1. The behavior of the function 
x{l) m this range is illustrated in Fig. 6. 

After expanding around the saddle point and performing the Gaussian integral over the 
small fluctuations, one obtains : 
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x(y) 



2.77 




0.5 



4.88 



y 



0.628 



Fig. 6. The BFKL eigenvalue x(l) an d the graphical solution to the saturation problem The 
value 7 S = 0.628 corresponds to the saturation saddle point and the value x{ls)/ls = 4.88 
determines the asymptotic energy dependence of the saturation momentum (for comparison, 
the "hard pomeron" saddle point 7p = 1/2 and its intercept cop = 41n2 = 2.77 are shown). The 
points p and A correspond to the proton and the nucleus respectively, for the same transverse 
momentum and rapidity, when they are far above the saturation scale. 



This is a reasonable approximation to Eq. (3.7) provided a s y ?Z 1, p ^ 1, and at least one 
of these quantities is much larger than one. 

When r varies from to oo, the function 70 (r) is monotonously increasing, and spans 
the range from 7o(0) = 1/2 to 70(00) = 1. Note, however, that physically the ratio r 
cannot be smaller than the value r s = p s (y)/a s y (with p s (y) = ln[Q;?(v4, y)/Q 2 s (A)}) which 
corresponds to the saturation momentum. Indeed, the linear approximation is justified 
only as long as k± > Q s (A,y), or p > p s {y). 

To compute the saturation momentum, and the corresponding saddle point 7 S = 7o(r s ), 
we shall use the saturation boundary condition (3.4), together with the approximation 
(3.12) for (pA(k±, y), which for this purpose is extrapolated down to k± ~ Q S (A, y). It can 
be justified that this is a legitimate procedure for the calculation of the energy dependence 
of Q s (A,y), although not necessarily also for its absolute normalization [38,39]. To the 
accuracy of interest here, Eqs. (3.12) and (3.4) imply just that the exponent in Eq. (3.12) 




(3.12) 
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must vanish on the saturation line: 



^(7o(r),r) 



7 s r s + x{ls) = 0, 



(3.13) 



r=r. 



which shows that r s is a pure number, to be denoted as c in what follows. This number 
together with 7 S are determined by Eqs. (3.11) and (3.13) as: 7 S = x{ls)/x'{ls) — 0.627... 
and c = t s = x\ls) = 4.883.... [7,37-39]. 

Note that in writing the saturation condition in the form (3.13) we have neglected the 
effect of the Gaussian fluctuations around the saddle point. The reason is that this effect 
is modified by the absorptive boundary condition at saturation [39] (that we neglect here), 
and thus cannot be correctly computed in the present approximations. In what follows, we 
shall often ignore such slowly varying prefactors, but only include the dominant behavior 
at high energy. To the same accuracy, the nuclear saturation momentum reads: 



This expression appears to satisfy the initial condition Q 2 (A,y = 0) = Q 2 (A), but this is 
only formal, since the calculation above is valid only for relatively large y, with a s y > 1, 
so the extrapolation to y — > is not justified. In fact, the energy dependence of Q 2 (A,y) 
displayed in Eq. (3.14) is only the first term in an asymptotic expansion at large y, out 
of which also the second term [39, 40] and the third one [41] are presently known. 

For practical purposes, the saddle point approximation (3.12) is still quite involved, since 
the function 7o(r) is known only implicitly, as the solution to Eq. (3.11). However, im- 
portant simplifications leading to fully explicit formulae can be performed either for very 
large momenta, p ^> d s y — this corresponds to the so-called "double-logarithmic accu- 
racy", or DLA, regime — , or for momenta above, but relatively close to the saturation 
scale, where one can expand Eq. (3.12) around the saturation saddle point. Below, we 
shall refer to the latter regime as the "BFKL regime" 6 . 

Let us briefly describe here the modifications which should be brought into the previous 
formulae in order to adapt them to the case of a proton. Recall that there is no saturation 
scale (i.e., no analog of the nuclear scale Q S (A)) in the proton at y — 0. The infrared 
cutoff is rather set by the scale Q p which is soft and reflects non-perturbative physics 
(cf. Eq. (2.18)). This is therefore also the scale which enters the formulae for the linear 
evolution of ip p (k±,y) with increasing y. E.g., the proton analog of Eq. (3.7) is obtained 

6 Note that this is not the usual BFKL regime associated with the saddle point tp = 1/2 and 
corresponding to the high energy limit a s y ^> p. Whereas 7p = 1/2 describes the evolution with 
y at fixed (and relatively low) k±, j s = 0.63 corresponds rather to an evolution where, when 
increasing y, kj_ is correspondingly increased, in such a way that the condition k\ > Q 2 (A,Y) 
(or p > p s {Y)) remains satisfied. 



Q 2 s (A,y) ~ Q 2 s (A)e' 



c = 



4.883.... 



(3.14) 
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by replacing Q S (A) — > Q p and taking an initial condition ^(7) = 1/(1 — 7) (compare to 
Eq. (3.8)). In particular, the variable p for the proton is defined as p(p, k±) = ln(/c^/Qp). 
Of course, for sufficiently large y, the gluon density in the proton will become large, and a 
(hard) saturation scale Q s (p, y) will be eventually generated. The dominant y-dependence 
of this scale at large y is described again by Eq. (3.14), but with Q S (A) — > Q p . More 
generally, all the formulae to be obtained below in this section can be translated to the 
case of a proton by replacing everywhere Q S (A) — > Q p . 



3.2.2 regime 

If p a s y, the saddle point is close to one, and only the first term (the one having a pole 
at 7 = 1) must be retained in the initial condition (3.8). A straightforward calculation 
yields 70 (r) ~ 1 — l/\/r, and : 



Q2 e V^~P 



Mk±, y) * — rr 1 , (3-15) 



where the factor Q\ has been reconstructed as Q 2 S {A)/ pa- This approximation is valid 
for a s yp 1, and as long as this condition is satisfied, Eq. (3.15) can be even used for 
relatively small y, with a s y <C 1. Indeed, in this regime Eq. (3.15) can be obtained from 
the DGLAP equation [57], which provides a resummation based on the assumptions that 
a s p > 1, but a s y <C 1. 



3.2.3 BFKL regime: weak universality and approximate scaling 

When p is larger than p s {y) = \n(Q 2 s (A, y)/Q 2 (A)) = ca s y, but not much larger, one can 
obtain an explicit form for <fA{k±,y) by expanding the general saddle point approximation, 
Eq. (3.12), in powers of r — r s (or, equivalently, in powers of (p — p s )/p s ). To second order 
in this expansion, one obtains [38-41] 



<p A (k ± ,y) ^^-{p- Ps + A} e-^~^ exp LiL-pl] , (3.16) 
a s N c { 2(3 a s y J 

or, more explicitly [recall that p — p s = \n{k\/Q 2 s (A,y))\ : 



Mk±,y) « {in yj^- + A) ( ^4^y S ex P 1-7^- 1^ T^K) (3-17) 

a s N c \ Q 2 (A,V) ) \ k 2 ± J F \ 2f3a s y Q 2 (A,v) J 
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In these equations, 7 S ~ 0.63 appears as an "anomalous dimension" (recall that, at very 
high k±, the actual power of l/k\ was one, both in the initial condition at y — 0, and in 
the DLA approximation (3.15) at y > 0). Furthermore, (5 = x"(ls) — 48.518... plays the 
role of a diffusion coefficient for the diffusion in \n(k\/Q 2 s ). 

The constant numbers k,q and A in Eq. (3.17) are not known 7 , as they depend upon the 
detailed matching onto the (unknown) solution in the transition region towards satura- 
tion. In particular, Eqs. (3.17) and (3.4) need not match with each other, since there is no 
reason why Eq. (3.17) should remain correct down to Q S (A, y). But in writing Eq. (3.17) 
we have nevertheless assumed that, when extrapolated to k± ~ Q s (A,y), the BFKL solu- 
tion becomes parametrically of the same order as the actual solution in the transition region 
(namely, they are both of order l/a s ). This may look surprising in view of the previous 
experience with the MV model where we have seen that, when extrapolated down to 
Qs(A), the twist piece ip^™ 1 remains parametrically smaller (by a factor pa) than the ac- 
tual distribution at saturation ip^. The two functions match with each other only at the 
(parametrically larger) momenta k\ ~ (In pa)Q 2 s {A). Still, as we shall see through a more 
detailed analysis in Sect. 5.3, the evolution with y is such that a power-law, 'twist', tail 
develops in the gluon distribution at momenta just above Q s (A,y), and this tail is pre- 
dominantly generated via radiation from those gluons which were originally at saturation 
(i.e., which at y = were included in tp*^). Thus, for a s y ^ 1, the twist contributions gen- 
erated by the evolution are not parametrically suppressed with respect to the saturating 
distribution anymore. 

An important property of Eq. (3.17), that we shall refer to as weak universality (as op- 
posed to the strong universality that has been found at saturation [35]) is the fact that 
the functional form of this expression, seen as a function of z = k\/Q 2 s (A, y) and y, is 
independent of the initial condition at y = (and, in particular, of A), but rather is fully 
determined by the quantum (BFKL) evolution together with the saturation boundary 
condition: <fA(k±,y) = <p(z,y), with ip an universal function. The initial condition affects 
only the value of the saturation momentum, cf. Eq. (3.14) [42]. 

The factor linear in p — p s in the right hand side of Eq. (3.16), or (3.17), does not come 
from the expansion of Eq. (3.12), but rather from the absorptive boundary condition at 
saturation [39,41], that we have ignored so far. The reason why this factor cannot be 
neglected (in contrast, e.g., with the prefactor in the expression (3.14) for the saturation 
scale) is that this factor introduces a dependence upon p which is not suppressed at large 
y. The other factors in Eq. (3.16), which involve p — p s in the exponent, do come from 
the expansion of Eq. (3.12). A priori, one would expect this expansion to be valid so long 
as (p — p s )/p s <C 1, since this is the "small parameter". However, since (3 is a relatively 
large number, the second order term in the expansion, quadratic in p — p s , remains small 

7 But it seems natural to impose A > 1/7, to avoid that Eq. (3.17) become a decreasing function 
of 1/Ar[ when approaching Q s (A,y) from the above. 
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as compared to the linear term even for p — p s ~ p s . We conclude that Eq. (3.16) (or 
(3.17)) is a good approximation for any p within the following range [38] : 

0<P-Ps(y)^p s (y), (3.18) 

often referred to as the "extended scaling window" , since within this range the BFKL so- 
lution approximately preserves the geometric scaling property characteristic of saturation 
[38, 39]: More precisely, Eq. (3.16) shows scaling so long as p — p s (y) <C \ / 2(3a s y, whereas 
for larger p, the scaling is violated by the diffusion term. Both the scaling property, and 
its violation, seem to be necessary to understand the small-x data at HERA [56,59,60]. 

Note finally that the existence of a relatively large window for geometric scaling above the 
saturation momentum is a consequence of the evolution 8 , so it takes some finite amount 
of rapidity (typically, a s y ~ 1) until this behavior is actually reached. 

3.2.4 The geometric scaling line : Approximately matching BFKL and DLA 

The approximate solutions that we have found so far apply in different kinematical regions 
- p ^> a s y for the DLA formula (3.15), and p— p s ^ p s for the BFKL expression (3.16) — , 
which have no overlap with each other. So, a priori, there seems to be no reason why these 
expressions could be matched onto each other. Still, if one chooses p = 2p s (y) = 2ca s y, 
and one estimates both the DLA and the BFKL expressions along this line, then one 
finds that the dominant (exponential) behavior in y is very similar, although not exactly 
the same, for both approximations. This is in agreement with an observation [38] that 
the geometric scaling line p = 2p s (y) is approximately the borderline between the regions 
dominated by the saturation saddle point and, respectively, by the DLA saddle point in 
the y — In k\ plane. 

This opens the possibility to achieve at least an approximate matching between the two 
limiting expressions, by appropriately tuning the borderline between these two regime. 
Of course, there is no fundamental reason in favor of this matching, but this will be 
convenient in practice, as it will allow us to avoid exponentially large discontinuities when 
studying the ratio 7Z p a in different regimes. Specifically, we shall tune the upper limit 
in the window (3.18) in such a way that the leading, exponential, behavior along that 
line be exactly the same for the DLA and BFKL approximations with fixed coupling. A 
straightforward calculation shows that Eq. (3.18) must be replaced by 

8 For instance, there is no such a window in the initial conditions at y = 0, where the property of 
geometric scaling can be defined with respect to the dependencies upon k± and A. As manifest on 
Eq. (2.13)-(2.15), it is only the saturating distribution which is just function of z = k\/Q 2 s (A); 
the twist terms have additional, explicit, dependencies upon A (and A), through factors of 1/pA- 
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< p-p s (y) < vp s (y) , 

or, equivalently, 



(3.19) 



Ql(A,y) « kl < Ql(A,y) = Ql(A,y)(^^y , (3.20) 
with v pa 1.708 determined by solving the following equation: 

/ (vc) 2 

c[(l - 7 > + 1] - V 4c (^ + 1) " = . (3.21) 

As anticipated, v is a number of order one. 

3.3 Linear evolution at high k± : Running coupling 

Since the (leading-order) BFKL equation corresponds to a fixed-coupling approximation, 
the inclusion of the running is a priori ambiguous. Physical arguments, together with 
recent experience [40,45] with the next-to-leading order BFKL equation [44], suggest 
that it should be appropriate to use a coupling which is running with the transverse 
momentum k± of the gluon [i.e., a s = a s (k\)}. This running has been already used 
in obtaining Eq. (3.2), valid in the saturation region. With such a running, the BFKL 
equation becomes difficult to solve in general, even formally, so it is interesting to keep in 
mind that, in order to study the high-energy limit 9 , one can also use the running with 
the saturation momentum [i.e., a s = a s (Q 2 s (A,y))}. Indeed, with increasing y, the relative 
separation from the saturation line (p — p s (y)) / Psiy) becomes smaller and smaller, and 
we can approximate, e.g. (with b = b N c /ir = 12N C / (11N C — 2N f )), 



s Mih^) ** ^ Q ' (A ^i 1 - w 1 ) • (3 - 22) 

In this equation r^(y) = ln(Q 2 (A, y)/A 2 ) = p s (y) + Pa (see Eq. (3.23) below for a more 
precise definition). By using just the first term in this expansion, i.e., d s (Q 2 (A,y)), one 
obtains the correct dominant behavior at high energy for both the saturation momen- 
tum and the gluon distribution [38-41]. However, interesting subleading effects at high 
energy (like the "diffusion term" in the gluon distribution) are not correctly given by 

9 In the present context, this is the limit in which, with increasing y, p is increased as well, in 
such a way to remain in the linear regime at p > p s (y)- 
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this approximation, whose applicability at finite y is restricted to a rather narrow strip 
(p — Ps) / Ps ^ 1 above the saturation line. To keep control on such subleading effects, we 
shall follow Refs. [39-41] and use the whole expression in the r.h.s. of Eq. (3.22) when 
studying the "BFKL regime" precursory of saturation. Of course, the resulting approxi- 
mation will be still limited to (p — p s )/p s <C 1 — because of the expansion performed in 
Eq. (3.22), and also of the other approximations needed when solving the (correspond- 
ing version of) BFKL equation [39,40] — , but as we shall see, the accessible window is 
considerably wider with this improved running. 

On the other hand, for much larger momenta, such that p — p s ^ p s , the running effects 
cannot be expanded anymore, and the actual running coupling a s (k\) has to be used. 
But in this regime, the BFKL equation itself boils down to the DLA equation, which is 
much simpler, and whose exact solution is known also with running coupling. 

As it should be clear from the previous considerations, the strategy for performing ap- 
proximations is now different from the fixed coupling case: Rather than solving first the 
general BFKL equation and then simplifying the result, it is more convenient to perform 
the approximations on the equation itself, and then solve a simpler equation. Since the 
approximations that we shall consider are different for the BFKL and the DLA regimes, 
our analysis will not cover the whole linear regime at p > p s (y) : the intermediate region 
between the BFKL and DLA regimes remains out of control. 

But before we present the general results, let us use the simplified running a s (Q 2 (A, y)) to 
give a rapid derivation of the leading energy dependence of the saturation momentum at 
large y. With this running, the only modification in the previous, fixed-coupling, analysis 
refers to the replacement d s y — > h A (y) in the exponent of the integrand in Eq. (3.5) (and 
the subsequent formulae). Here, h-A^y) is defined as the solution to dtiA/dy = a s (Q"j,(A, y)) 
with /ia(0) = 0. Then calculations entirely similar to those described in the previous 
subsection lead to the condition r s = p s {y)/h.A{y) = c, where c = 4.883... and p s (y) = 
\n(Q 2 (A, y)/Q A ), as before. This condition together with the equation defining tiA(y) 
form a coupled system of equations for the functions tiA^y) and Q 2 (A,y), whose solution 
is conveniently written as [38,39,42]: 



For completeness, we note that ch A (y) = ta(u) — Pa- 

For relatively small y, such that 2cby <C p 2 A , one can expand r^(y) — pa + ipcy/pA)- Then, 
Eq. (3.23) reduces to the same expression as for fixed coupling, cf. Eq. (3.14), but with 
a s evaluated at the initial saturation scale (i.e., a s = a s (Q A ) = b/ pa) [40]. 





(3.23) 
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But with increasing y, the dependence on A becomes weaker and weaker 10 [42] 



Q 2 s (A,y) ~ A 2 ev^ exp |^^} for 2cby » . (3.24) 

By using similar formulae (with pa — > p P ~ 1) also for the proton, it becomes clear 
that, unlike in the fixed coupling case, where the ratio Ql(A, y)/Q 2 s (p, y) = Q 2 S (A) / Q 2 s (p) 
is independent of y and thus fixed by the initial conditions, in the case of a running 
coupling, this ratio decreases with y, and converges to one at large y [42] : 



i£fH- {M}- 1 for *»* A - (x25) 

Thus, due to running coupling effects in the evolution, the initial difference between the 
proton and the nucleus is washed out after a rapidity evolution cby ~ p\. 



3.3.1 DLA regime 

With a running coupling a s {k±), the solution to the DLA equation — the simplified 
equation valid at high kj_, which is a common limit of the BFKL and DGLAP equations 
- is well known, even exactly. For consistency with the other approximations that we 
have considered so far, we shall only use the corresponding saddle point approximation, 
valid when 4byr] ^> 1, with r\ ~ hip (see below for a precise definition). This reads: 



<PA(k±,y) ~ t\ttt 



' Abyr] 



In 



In 

111 A 2 
111 A 2 



We expect this to be a good approximation so long as p - 
practice we shall use under the stronger form p — p s > t~a- 



(3.26) 

p s > p s , a condition that in 



3.3.2 BFKL regime 

The solution to the BFKL equation with the running coupling expanded as in Eq. (3.22) 
has been obtained in [39] (see also Refs. [40,41]) for the case of the dipole scattering am- 
plitude. In what follows we shall adapt this solution to the case of the gluon distribution, 

10 Choosing pa = 4 (a reasonable value for the RHIC phenomenology) and using 2cby ~ lOy, 
one finds that the running effects become important — in the sense that the two terms under 
the square root in Eq. (3.23) become comparable — already for y of order one. 
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explain its properties, and carefully establish its validity limits (thus going beyond the 
original discussion in Refs. [39-41])). 

When applied to the gluon occupation factor, the solution in Refs. [39-41]) reads 

« ■ ^^<^^H-^)' (3 - 27) 

with the following notations: Kq is a pure constant of 0(1), ta{d) has been defined in 
Eq. (3.23), 7 S = 0.628 is the saturation saddle point, z = k\/Q 2 s (A,y), and therefore 
hiz = ln(kj_/Q 2 s ) = p — p s . Furthermore, Ai(x) is the Airy function, £i = —2.33... is its 
first (rightmost) zero, A is a constant that we are not able to determine, 



D=[x''( 7s )/(2x( 7s ))] 1/3 = 1.99..., (3.28) 

is the "anomalous" diffusion coefficient, and (3 = x"(ls) — 48.5... is the "standard" one 
(however, notice the 3/2 factor in the Gaussian in Eq. (3.27), instead of a 2 in Eq. (3.16)). 
Note that, in the present context, we find it more natural to measure the transverse 
momenta with respect to A, rather than with respect to the original saturation scale (as 
we did for fixed coupling). 

To the same approximation, the saturation momentum is given by 

Ql(A, y) = A 2 exp [r A - ^ rf + a) . (3.29) 

Since obtained after high energy approximations, Eqs. (3.27) and (3.29) apply only for y 
large enough, such that by/pA > 1- (This is the usual condition d s y > 1 for the validity of 
the BFKL approximation, where a s is now evaluated at the original saturation momentum 
Qa-) As expected, Eqs. (3.29) and (3.23) show the same dominant behavior at large y. 

The validity range of Eq. (3.27) in z is determined by the various approximations per- 
formed in deriving this result [39]: First, there is the condition Inz <C ta, which has two 
sources: (i) the expansion of the BFKL x function around the saturation saddle-point 11 
(the analog of the saddle-point approximation in the fixed coupling case), and (ii) the ex- 
pansion of the running coupling around the saturation line, cf. Eq. (3.22). Second, there 
is a more restrictive condition, related to the further approximations performed when 
solving the "diffusion equation" (Eq. (74) in Ref. [39]) obtained after expanding the x 



11 By itself, this expansion could be justified up to Inz ~ ta, because of the relative large value 
of the coefficient [3. 
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function. This condition leads to the following window for the applicability of Eq. (3.27) : 

< lnz = p-p s < ^5- rf , (3.30) 

V5 

to be referred as the 'BFKL regime' in the presence of a running coupling. This window 
includes several interesting kinematical regimes, as we discuss now: 

Using the expansion of the Airy function around its zero, 

Ai(fi + a;) =Ai'(fi) (^ + J x3 + '- ^ for a; < 1, (3.31) 
one can deduce an approximation for ip valid near to the saturation line, namely, 

Mk±,y) « - z~ ls (Inz + A) , (3.32) 

a s (k±)J\ c 

where trivial constants have been absorbed in the unknown multiplicative factor k,q. Ex- 
cept for the mild scaling violations brought by the factor l/a s (k±), this expression is a 
scaling function, i.e., a function of z — k\/Q 2 s (A, y). Eq.(3.32) is valid so long as the cubic 
term in (3.31) is smaller than the linear one; this condition implies k± <C Q g (A,y) where: 



Q 2 M,y) » Ql(A,y) exp (J^^T) . (3.33) 

Thus, it is natural to choose Q g (A,y) as the upper boundary for the region of geometric 
scaling in the present case. Note that, with running coupling, this region is only a small 
part of the BFKL regime defined by Eq. (3.30). Since ta(v) is monotonously increasing 
with A and y, and approaches an A-independent value when y — > oo, it is clear that 

Q 2 g (A 1 ,y)>Q 2 g (A 2 ,y) for A 1 > A 2 , (3.34) 

and the two scales will approach each other at very large y. Note also that the expression 
in Eq.(3.32) is a universal function, in the sense of Sect. 3.2.3. 

Consider also the region outside the scaling domain, at Inz > r^ 3 . Using the asymptotic 
form of the Airy function, 

Ai(x) w ^ exp [~ x 3 ^ for x > 1, (3.35) 

one immediately finds 

. . tT 2 ( , 2 ln 3 / 2 z 2 \n 2 z\ , , 

w ^^h7^ ex H" 7s ^"^"^^J' (3 - 36) 
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where Ko has been redefined in order to absorb some trivial constants, and A has been 
dropped since it is small compared to In z in the region of interest. Note the appearance 
of diffusion terms of two types — from the Gaussian Eq. (3.27), and from the asymptotic 
expansion of the Airy function — , which both lead to scaling violations. There is also a 
mild violation from the prefactors. 

Moreover, at variance with the fixed coupling case (cf. Eq. (3.16)), the scaling violations 
in Eq. (3.36) are synonymous of universality violations : indeed, the "diffusion time" ta(v) 
(which replaces the a s y of the fixed coupling case) is now explicitly dependent on A. But 
for sufficiently large y, the A-dependence of ta(v) becomes very weak, as we have seen, and 
universality gets restored when 2cby 3> p\. For even larger y, such that 2cby ^> p\, the 
saturation momentum becomes independent of A, cf. Eq. (3.25), so the gluon occupation 
factor — which in this regime is a universal function of z and y — is then the same in 
the proton and in the nucleus for any k± in the BFKL regime. 



4 High-A; ± suppression from quantum evolution: The general argument 

Before we embark ourselves in a systematic analysis of the ratio 1Z p A(k±,y), let us present 
here a general argument explaining why the quantum evolution leads to a suppression in 
this ratio at generic momenta. The argument can be succinctly formulated as follows 12 : 
The suppression of the ratio 1Z p A(k±,y) with increasing y is due to the different evolution rates 
for the gluon distribution in the nucleus and in the proton: the proton distribution evolves faster 
because (a) the corresponding saturation momentum in the initial conditions is much smaller 
(Q p <C Q S (A)), and (b) the convexity of the function x(l) (the eigenvalue of the BFKL 
kernel) accelerates the evolution with increasing p. In particular, the difference between the 
proton and the nuclear evolution rates is particularly pronounced at small y (when the 
proton is in the DLA regime for the interesting values of k±), which explains the rapid 
suppression seen in the early stages of the evolution [25]. 

Moreover, similar arguments will allow us to conclude that, after some initial evolution 
a s y ^ 1, and for transverse momenta in the proton perturbative region (k± ^> Q s (p,y)), 
the ratio TZ p A(k±,y) is monotonously increasing with k±, and a decreasing function of A. 
The monotonic behavior with k± shows that the Cronin peak has already disappeared 
for such values of y. As for the A-dependence, this can be related to the experimental 
observation that the central-to-peripheral ratio Rep for (pseudo)rapidity rj > 1 is smaller 
for central collisions than for more peripheral ones [4-6]. (Peripheral collisions involve less 
nucleons than the central ones, thus effectively corresponding to smaller values of A.) 
This tendency is opposite to that identified at y = in our theoretical calculations (cf. 

12 We would like to thank Al Mueller for helping us clarifying this general argument. 
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Eqs. (2.23)-(2.24)), and also at r\ = in the experimental results at RHIC [4-6], where 
the amplitude of the Cronin peak appears to be increasing with A. 



In this context, it is interesting to notice that the quantity a (the logarithmic slope in 
the parametrization Rqp ~ e av ) which has been measured at BRAHMS [4] corresponds 
in our approach to the right hand side of either Eq. (4.6), or Eq. (4.13), for the nucleus 
in the linear regime, or in the saturation regime, respectively. With this interpretation, a 
is dominated by the linear evolution of the proton, and is not the same as the saturation 
exponent (in contrast to a suggestion in Ref. [4]). 

The general arguments to be presented in this section apply for the case of fixed coupling 
alone: for a running coupling, we have used different approximations in the BFKL and 
DLA regimes, so, based on these approximations, one cannot construct a general argument 
valid in the whole perturbative region. But in Sect. 6 we shall check through explicit 
calculations that all the conclusions about the general behavior of the ratio 1Z p A(k±,y) to 
be established below in this section remain true for a running coupling. 

4-1 A general argument on the y -dependence 

The ^/-dependence of the ratio in Eq. (1.1) is encoded in: 

dlnTZpA dlnipA d\nip p 
dy dy dy 

We shall study separately the linear and the saturation regimes for the nucleus. The 
proton will be always in the linear regime, as described by the BFKL solution in the 
saddle point approximation, Eq. (3.12). 

4-1.1 Nucleus in the linear regime (k± ^> Q s (A,y)) 

In this regime, both the proton and the nucleus are described by Eq. (3.12), which implies: 



where we have neglected the slowly varying prefactors. Here, we have introduced the more 
specific notation: 



ln(p A (k ± ,y) 



a s yF('jo(r A ),r A ) , 



(4.2) 



p(A,k ± ) _ \nk 2 JQ 2 s (A) 



(4.3) 



oi s y 



a s y 
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(with Q S {A) — > Q p for the proton) to emphasize that quantities like p = p(A, k±) and ta 
are different for the proton and the nucleus, even when considered for the same values of 
the kinematical variables k± and y. Since Q S (A) 3> Q p then, clearly, 

p(A, k ± ) < p{p, k±), and r A (k ± ,y) < r p (k±,y) . (4.4) 

After taking a total derivative w.r.t. y in Eq. (4.2), and using the saddle point condition, 
Eq. (3.11), together with the definition (3.10) of the function F, one finds : 

—, ~ u s x{ia) , (4.5) 

dy 

(with the simplified notation 7^ = 70 (r^)) and therefore: 

ldlnn ? A « X (7A) - X(7 P ) , (4-6) 



«5 C?l/ 



which is always negative (for and y within its range of validity), as we show now: Since 
x{l) is a convex function (x"(t) > for < 7 < 1), the function 70 (r) which gives the 
saddle point is monotonously increasing: 

= I > (4 7) 

dr x"(7o(r)) >U ' 10 

(this follows by differentiating the saddle point condition (3.11) w.r.t. r), which together 
with Eq. (4.4) implies that 7^ < 7 P , from which Eq. (4.6) finally follows because x(t) 
is monotonously increasing for any 7 > 1/2 (recall that 7^ > 7 S ~ 0.63). Some of these 
properties are manifest on Fig. 6. 

Moreover, the suppression rate (4.6) is largest in the early stages of the evolution, but 
decreases with y, and approaches zero (through negative values) when y — > 00. Indeed, 
the difference: 

olsV ot s y a s y 

is independent of k± and largest at small y, but it decreases with y, and vanishes asymp- 
totically at very large y. Therefore, when y — > 00, we have 7^ — 7 P — * 0, and thus also 
x(1a) ~ Xilp) -»■ 0. 

To be more specific, consider two interesting physical regimes: 
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a) For very large k± ^> Q g (A, y), both the proton and the nucleus are in the DLA regime, 
where r >> 1, 7o(r) ~ 1 — x(7o) ~ \A> an d therefore: 



1 ciln^pA _ lp(Ak±) I p{p, k±) ^ ^ 



a s dy V V a sV 

which confirms that the rate of variation is largest in the early stages of the evolution 
(although Eq. (4.9) cannot be trusted, strictly speaking, when y — > 0). 

b) For kj_ and y such that both the nucleus and the proton are in the BFKL regime 
(which becomes possible when y is large enough for Q g (p, y) > Q S (A, y) ; see Fig. 1), then 
one can expand both x(1a) and x(lp) around 7 S fa 0.63. This is entirely similar to the 
expansion leading from Eq. (3.12) to Eq. (3.16), and gives: 



1 d\nK pA „ c PA / 1+ Wjfcj + Pa\\ (41Q) 



a s dy I3a s y { ca s y \ Q 2 s (A,y) 2 

which shows that, when increasing y and k± in such a way that \a{k\/ Q 2 S {A, y)) ~ const., 
the rate eventually becomes independent of k±, and vanishes as 1/y. 



4-1.2 Nucleus in the saturation regime (k± < Q s (A,y)) 

All we need to know about this regime for the present purposes is that the nuclear gluon 
distribution at saturation is only a function of z = k\/Q 2 s {A, y). To remain in the satura- 
tion region when increasing y, one must simultaneously increase k± along a line parallel 
to the saturation line, such that z remains fixed. Proceeding as before, we have: 



zQ 2 (A y) 

lnn pA {z,y) ^ lmp A (z) - a s yx(l P ) +7 P hi — s — 1 — , (4.11) 
where the proton saddle point 7 P is now determined by (cf. Eq. (3.11)) : 



\n[zQ 2 (A,y)/Q 2 P } HzQ 2 (A)/Q 2 P } , , . , ^ 

X (7p) = ^Ty = ^y + X (7s) - (412) 

In the above equation, the second equality follows after using Eq. (3.14) for the saturation 
momentum and recalling that c = x'ils) ~ 4.88 with 7 S fa 0.63. Taking a derivative w.r.t. 
y in Eq. (4.11) and using Eq. (4.12) yields: 
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1 dhilZpA 



a s dy 



~ ipx'ds) - xiip) , 

z 



(4.13) 



which is negative for any y, as is quite obvious by inspection of Fig. 6 : The straightline 
f(j p ) = jpx'ils) is the line which is tangent to x{lp) a ^ lp — Is] but x(lp) is a convex 
function, and thus it lies above its tangent for any 7 P ^ j s . Moreover, the suppression 
rate is very large at small y (since xilp) ~~ > 00 as lp ~^ 1)> but it decreases to zero when 
y — > oo (recall that 7sX'(7s) — x(7s))- We conclude that the ratio lZ p A(z,y) is rapidly 
decreasing with y in the early stages of the evolution, and approaches a constant value 
for very large y. 



4-2 A general argument on the k±- dependence 



We now present a mathematically similar argument which allows us to characterize the 
k± dependence of the ratio 7Z p A(k±,y). We would like to show that, for a fixed rapidity 
y with a s y ^ 1 (such that the saddle point approximation (3.12) applies), and for all 
momenta k± ^> Q s {p,y) (where the proton is in a linear regime), the function TZ p A(k±,y) 
is monotonously increasing with k±. 

Consider first the linear regime for the nucleus, at kj_ ^> Q S (A, y). Using Eq. (4.2) and 
the saddle point condition, Eq. (3.11), one obtains (e.g., for the proton), 



dln ^ „ _ 7 U 14) 

dln(fci/A 2 ) 7p ' 1 } 

together with a similar formula for the nucleus, and therefore: 



dlnKpA 1K v 

which by the previous arguments is strictly positive and vanishes asymptotically when 
k±_ — > oo. Thus, in this regime, the function TZ p A(k±,y) is strictly increasing with k± for 
any fixed y, and saturates to a constant value when k± — > oo. As it will be checked in 
Sect. 6, this limiting value is one. 

Consider also the case where the nucleus is deeply at saturation (k± Q s (A,y)). Then, 
ifA is only logarithmically decreasing with z = k\/Q 2 s {A, y) (see Eq. (3.1)), while ip p has a 
faster, power-like, decrease, with "anomalous dimension" 7 P > 7 S ~ 0.63 (cf. Eq. (4.14)). 
Therefore, in this range too, the function 1Z p A(k±,y) is increasing with k±. 
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The previous analysis does not cover the nuclear transition region at k± ~ Q S {A, y). But 
the discussion following Eq. (3.17) shows that, when a s y ?t 1, there is no large (i.e., 
parametrically enhanced) mismatch between the value of ipA at k± ~ Q s (A,y), Eq. (3.4), 
and the corresponding extrapolation of the BFKL solution (3.17): the original mismatch 
between cp^ 1 and </?^ lst at k± ~ Q S (A) (cf. Sect. 2.2) has been washed out by the evolution. 
This property, together with the above finding that ipA is monotonously increasing with k± 
on both sides of the saturation line, precludes the existence of a (parametrically enhanced) 
peak at k± ~ Q s (A,y). That is, the peak must have already disappeared when a s y 1. 
The mechanism for this disappearance will be clarified in Sect. 5.3. 



4-3 A general argument on the A-dependence 



Finally, let us use similar arguments to study the dependence of the ratio 1Z p a upon the 
atomic number A. We shall find that, after an evolution 13 a s y ^ 1, the ratio becomes 
a decreasing function of A. We shall consider the A-dependence for two cases: (a) with 
fixed k± and y (this corresponds to the real experimental situation for the measurement 
of R C p at RHIC) and (b) with fixed z = k\/Q 2 s (A,y) and y (as relevant, e.g., for the 
A-dependence of the Cronin peak). 

Notice first that the dominant A-dependence enters the gluon occupation factor through 
the ratio between k\ and the saturation momentum (see, e.g., Eqs. (3.1) and (3.12); the 
factor Xp~A in the latter brings only a subleading, logarithmic, dependence on A, whose 
contribution to the derivative in Eq. (4.16) below is suppressed by a factor 1/pa)- 

(a) With fixed k± and y, <p p (k±,y) is clearly independent of A, and one can rely on 
Eq. (3.12) (nucleus in the linear regime) to write: 



d\mp A 9 In A 1 / 3 dp din ip a d In A 1 / 3 
. dA dA = dA~dp dA 

-{-ia) - -7TT = -^rri 1 - 7a) < 0, (4.16) 



dA y ' ' dA dA 



where we have also used p(A, k±) = \nk\/ Q 2 S (A) = In k\/A 2 — pa and the relation (4.14) 
written for the nucleus. The same conclusion holds when the nucleus is deeply at satu- 
ration, since there ipA is only logarithmically varying with kj_/Q 2 (A,y) (and thus with 
A). Finally, for k± around Q s (A,y) one can rely on the same matching argument as in 



13 This is a conservative estimate, necessary for the validity of the general arguments in this 
section. But the explicit calculations in Sect. 5.1 will reveal that, for kj_ ~ Q s (A,y) at least, 
1Z P A starts to be a decreasing function of A already after the shorter evolution a s y ~ 1/pa- 
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the previous discussion of the /c_i_-dependence to conclude that the ratio 1Z p a preserves a 
monotonic, decreasing, behavior with A. 



(b) Consider now the situation in which A and k± are increased simultaneously, in such 
a way that z = k\/Q 2 s (A,y) stays fixed. Then p(A,k±) = \nk\/Q 2 s {A) is fixed as well, 
and the gluon distributions of the nucleus (seen now as a function of z and y) has no 
extra dependence on A. On the other hand, the distribution of the proton becomes now 
dependent on A, through k±: Writing the proton p- variable as 



pip, k±) = In 



i,2 



In z + In 



~ Inz + ca s y + p A , 



(4.17) 



Ql Ql 

where the last term carries the dependence on A, and using Eq. (4.14), one deduces that 



d In ip p 



dA 



-1p 



z,y 



dp A 
OA 



(4.18) 



This yields the expected conclusion that 1Z p a is decreasing with A : 

dp a, 



M ln7ZpA 



z,y 



OA 



(1-7p)<0. 



(4.19) 



Note also that the rate of decrease in 1Z p a is faster when measured with k± and y fixed 
than with z and y fixed (because 1 — 7^ > 1 — j p ). 



5 The evolution of the Cronin peak with increasing y 



The problem of the evolution of the Cronin peak with increasing y is a delicate one, 
as it involves the nuclear transition region towards saturation, on which we have little 
analytic control at y > 0. Moreover, the experience with the MV model suggests that 
the existence of the peak and its properties are tributary to the actual behavior of the 
nuclear gluon distribution tpA{k±,y) at momenta just above the saturation scale: If the 
distribution is rapidly decreasing with k± (say, according to an exponential law), then a 
pronounced peak exists; but if its decrease is only power-like, then there is at most a very 
flat peak, if any. In the MV model, the gluon distribution above Q S (A) has been found to 
be the superposition between an exponential and a power-law tail, with the exponential 
being the dominant contribution though, since parametrically enhanced at large A. As 
we shall see below in this section, the effect of the evolution is to enhance the power-law 
contributions, which for a s y ~ 1 become as large as the exponential one (for momenta 
just above Q S (A)). When this happens, the Cronin peak has completely flattened out. 
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In order to study this flattening, one needs a more accurate calculation which follows 
the non-linear evolution of the nucleus in the saturation region. To achieve the necessary 
accuracy while still preserving an analytic control, we shall perform just one step in the 
evolution of the MV model according to the non-linear Kovchegov equation [33]. Strictly 
speaking, this calculation applies only for rapidities a s y <C 1, but as we shall see in Sect. 
5.3, this is enough to reveal the mechanism responsible for the flattening of the peak. By 
extrapolation, we shall then conclude that the peak disappears after an evolution a s y ~ 1. 

Whereas the flattening of the peak is a rather subtle effect which has to do with the 
nuclear evolution, the rapid decrease in the height of the peak, on the other hand, is 
a more robust phenomenon, and also easier to calculate, since this is due solely to the 
perturbative evolution of the proton. Because of that, we shall begin our analysis in this 
section with a study of the evolution of the magnitude of the peak with increasing y. 
Since we expect the peak to follow the nuclear saturation momentum, we shall consider 
the evolution of the ratio 1Z P A(k±,y) along the nuclear saturation line k± = Q s (A,y). 
This is possible within the present approximation because, along this particular line, the 
nuclear gluon distribution is known, cf. Eq. (3.4). By itself, this calculation cannot tell us 
whether a peak actually exists or not; but so long as the peak exists, it provides us with 
a correct estimate for the magnitude of the peak and its parametric dependencies upon 
A and y. Although, as we shall argue later, the peak disappears already when a s y ~ 1, 
in Sects. 5.1 and 5.2 we shall follow the evolution along the nuclear saturation line up to 
very large y. This will reveal the basic features of the high-/cj_ suppression, to be more 
systematically analyzed in Sect. 6. 



5. 1 The suppression of the peak: Fixed coupling 



In this subsection and the following one, we shall need only the (proton and nuclear) 
gluon distributions along the nuclear saturation line k± = Q s (A,y). For the nucleus, 
this is simply a constant times l/a s , cf. Eq. (3.4). For the proton, the corresponding 
distribution is in a linear regime, which can be either DLA, or BFKL, depending upon 
the value of y: With increasing y, the proton geometric scale Q g (p, y c ) rises faster than the 
nuclear saturation momentum, so the corresponding evolution lines in the plane y — In k\ 
cross each other at some rapidity y c (see Fig. 1), where the proton changes from the DLA 
to the BFKL regime. The condition Q s (A,y c ) = Q g {p,y c ) together with Eqs. (3.14) and 
(3.20) imply : 



1 Q S {A) p A 

ca s y c ~ - In *V ~ — . 5.1 
v Q 2 v 
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5.1.1 y < y c : Proton in the DLA regime 



For y < y c , but such that a s yp > 1, one can use Eq. (3.15) with Q S {A) — > Q p to deduce : 

K sat (A,y) = K P A(k ± = ft(A2/),2/)~f^e-« 

~ p A exp|ca s y - ^4a s y(p A + ca s y) J, (5.2) 
where we have used Eqs. (2.7), (2.19) and (3.14) to write: 

ff(^y) _ QiiAy) _ asy Q 2 S {A) _ 
and also (recall that Q 2 V ~ A 2 ) : 



p ee In = ca s y + In ~ c « sy + pA . (5.4) 

Note that, in writing Eq. (5.2), we have ignored the slowly varying prefactor in Eq. (3.15), 
since presently we are only interested in the dominant parametric dependencies upon A 
and y. When y — ■> 0, Eq. (5.2) is formally consistent with the corresponding result in the 
MV model, Eq. (2.22), but one should keep in mind that the results obtained here cannot 
be used for very small values of y. 

Eq. (5.2) offers, in particular, an estimate for the magnitude of the Cronin peak; indeed, 
for as long as it survives, the peak should be located in the vicinity of Q S (A, y). Note that, 
for any < y < y c , the expression within the exponent in Eq. (5.2) is negative, showing 
that the magnitude of the peak is rapidly decreasing with y. Since the initial maximum 
at y = was relatively large (1Z max (A) ~ p A , cf. Eq. (2.23)), it is interesting to check how 
fast is the height of the peak decreasing to a value which is parametrically of order one. 
The condition lZ sa t(A, y) ~ 1 for y ~ y implies: 



ln 2 p A (lnlnA 1 ^) 2 
a s y - ^- « a s y c . (5.5) 

This is a small rapidity interval, but is still within the reach of the saddle point approx- 
imation (3.15), since pA(x s yo ~ In 2 Pa > 1- In fact, for such a short evolution in y, one 
can ignore the evolution of the nucleus in the vicinity of the peak, i.e., one can neglect 



55 



ca s y <C 1 next to pa within Eq. (5.2). This is reassuring since, when a s y <C 1, one cannot 
really trust the BFKL approximation (3.14) for the saturation scale of the nucleus. 

In view of this, it is possible to obtain a more accurate estimate for the rapidity y by 
using the MV model for the gluon distribution of the nucleus together with the complete 
DLA expression (3.15) for the proton distribution, including the prefactor. This yields: 



a = V2^z T(0, zq) ~ 0.706 . (5.6) 

The corrections to this result (coming either from the evolution of the nucleus, cf. Sect. 
5.3, or from the twist terms in the MV model, that we have neglected here) are of 0{p~ A 2 ). 

Eq. (5.2) also shows that, after a short rapidity evolution a s y ~ 1/pa, the magnitude of 
the Cronin peak becomes a decreasing function of pa (and thus of A), in sharp contrast 
with the corresponding behavior at y — 0, cf. Eq. (2.23). This behavior is in agreement 
with the general arguments about A-dependence given in Sect. 4. 

As anticipated, the rapid decrease in the height of the Cronin peak during the early stages 
of the evolution is to be attributed solely to the fast evolution of the proton. For larger 
rapidities, y 3> yo, the evolution of the nucleus starts to matter as well (and Eq. (3.14) 
can be trusted), but as long as y < y c the proton evolution is still faster, so lZ sa t keeps 
decreasing with y, as manifest on Eq. (5.2). In particular, for y = y c one has: 

K sat (A,y c ) ~ PA e- KpA ~ « 1, k = \ ( ^4(1 + u)/c - l) » 0.29 . (5.7) 

For even larger y, Eq. (5.2) would predict an increase of the peak with y, but this is 
incorrect, since for y > y c this equation does not apply any longer. 

5.1.2 y > y c : Proton in the BFKL regime 

For y > y c , the proton enters the scaling window (3.20), where Eq. (3.17) becomes ap- 
propriate. We shall shortly argue that, when this happens, the Cronin peak has already 
disappeared; but it is still interesting to follow the ratio 1Z P A(k±,y) further up along the 
nuclear saturation line. By using Eq. (3.17) with Q S (A, y) — > Q s (p,y), together with the 
relations: 

= 2^ = ^, M-4V3 P ^, 4 , ( , 8) 



a s yo 



4:PA 



In 2 



ap A Jln(ap A ) 
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one finds (with the simpler notation 7 = j s ) : 



Note the emergence of the power of A in the denominator, which provides a strong sup- 
pression factor which is independent of y. This power was missing in the DLA (compare 
to Eq. (5.2)), but appears here as a consequence of the 'anomalous dimension' 7 < 1 
characteristic of the BFKL solution in the vicinity of the saturation line [38,39]. 

Note furthermore that Eq. (5.9) shows only a weak dependence on y, coming from the "dif- 
fusion term" in Eq. (3.17): the dominant dependencies, which are exponential, have can- 
celled in the ratio between the nuclear and the proton saturation momenta, cf. Eq. (5.8). 
Thus, as compared to the DLA regime at y < y c , where lZ sa t(A, y) is rapidly decreasing 
with y, in the BFKL regime at y > y c this ratio almost stabilizes at a very small value, 
proportional to an inverse power of A. 

In particular, when approaching y c from the above, one finds: 



W V ' ] ~ (I^f ~ (AUHnA^f (5 ' 10) 

where k has now been generated as k = 1 — 7 — vc/2(3, but this is the same power 
as in Eq. (5.7), by construction (because of our matching between the DLA and BFKL 
approximations along the geometric scaling line). The subleading dependencies on A, 
which are logarithmic, are different in Eqs. (5.7) and (5.10), but this difference merely 
reflects our imperfect matching. 

Furthermore, for very large y, such that 2(5a s y ^> p\ (in practice, this limit is approached 
very fast, since 2f3 ~ 97 is a large number), the ratio becomes independent of y : 



Kat(Ay) - ( A i/3 pA )i- 7 ~ (AV3i n AV3)i-7 for » Pa- (5.11) 



To summarize, the inverse power of A which characterizes the suppression increases from 
k ~ 0.29 at the end of the DLA regime to 1 — 7 ~ 0.37 at asymptotically large y, which 
shows that a large fraction of the suppression is actually built during the DLA phase of 
the evolution. 
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5.2 The suppression of the peak: Running coupling 

With running coupling, the kinematical domains in which our 'DLA' and 'BFKL' approx- 
imations apply are well separated from each other, and the intermediate region between 
the two regimes is not under control. From Sect. 3.3, we recall that for very large k±, such 
that Inz > T p (y) (with z = k\/Q 2 s (p,y)), the proton is described by the DLA approxi- 
mation (3.26), while for intermediate k±, within the range < Inz < const x r 2 / 3 (y) (cf. 
Eq. (3.30)), the BFKL approximation (3.27) applies. Thus, when the momentum is taken 
along the nuclear saturation line (k± = Q s (A,y), or Inz ~ ta(v) — T P (y)), the proton is 
in the DLA regime for y < y c , and in the BFKL regime for y ^ y' c . Here, 2cby c ~ p 2 A , and 
y' c is determined from the condition Ta(v) — T P (y) ~ const x r 2 ^ 3 (y) as: 



where the constant factor in the r.h.s. has been obtained according to Eq. (3.30). For 
even larger y > y", the proton enters its geometrical scaling region, where the simpler 
expression (3.32) applies; y" is clearly determined by the condition Q 2 (A,y") = Q 2 (p,y"), 
which together with Eq. (3.33) implies: 



Finally, for y so large that 2cby ^ p\, the proton saturation momentum catches up with 
that of nucleus, cf. Eq. (3.25), and then the proton too enters a non-linear regime. The 
constant factors appearing in the above equations for y c , y' c and y" c are not really under 
control, but their parametric dependencies on p A are correct as written. 

At this level, several differences with respect to the fixed coupling case are already mani- 
fest: First, there is an interval in rapidity that remains out of control, namely 



Second, the characteristic rapidities y' c and y", at which the proton enters the BFKL and 
scaling regimes respectively, are parametrically larger than the (unique) corresponding 
scale for fixed coupling, namely y c in Eq. (5.1). This shows that the evolution is slower in 
the presence of a running coupling. Third, with a running coupling, the linear (BFKL) 
approximation for the proton does not apply up to arbitrarily large y, but only so long 
as 2cby < p\. But for larger y, and k± = Q s (A,y), the proton and nuclear distributions 
coincide with each other, so the ratio 1Z sat reduces to the trivial factor 1/A 1 ^ 3 . 

Other differences between the fixed and running coupling scenarios will be revealed by 
the subsequent analysis. 




(5.12) 




(5.13) 



Pa Z 2c% < p X l 



(5.14) 



58 



5.2.1 y < y c : Proton in the DLA regime 

In this regime, the proton is described by Eq. (3.26) where, for k\ = Q 2 (A, y), the variable 
rj can be evaluated as (see also Eq. (3.23)): 

rj = \ n \ n 9li^yl _ lnln?| ~ hi r A (y) . (5.15) 
A 2 A 2 

By also using Eq. (3.4) for the nucleus, one finds : 

TZ sat (A,y) ~ r A (y) exp|r A (y) - p A - \j 4by\xiT A (y)^ (5.16) 

which for y — > is of 0(p A ), but it decreases rapidly with y, and it becomes of order one 
already after a small rapidity evolution y , with: 

y ~ 1 lnp A ~ In A 1 / 3 . (5.17) 

The saddle point approximation for ip p (k±,y), Eq. (3.26), is justified for y ~ y Q since 
^4by r] ~ ^/Wy^hTp~ A = In p A > 1. We find a situation similar to the fixed coupling case 
(cf. Eq. (5.5)), namely, yo is so small that the corresponding evolution of the nucleus can 
be neglected: r A (y ) — p A ~ In p A j p A 1, and therefore Q 2 S (A, y) ~ Q\ for y < y . In 
particular, one can use the MV model for the nucleus in order to obtain a more accurate 
estimate for y , valid up to corrections of 0(p A 2 ) (the analog of Eq. (5.6)) : 



y ° = TkT — ln2 

46 ln p A 



ap A J\n(ap A ) 



a = v^oI^O, zq) ~ 0.706 . (5.18) 



When comparing Eqs. (5.5) and (5.17), it appears that the running coupling estimate for 
yo is parametrically larger than the fixed coupling one in the limit where A is large. This 
confirms the fact that the running effects slow down the evolution. However, one should 
keep in mind that the rapidity is naturally measured in units of l/a s , and whereas a s yo ~ 
In 2 p A / p A for fixed coupling, we also have a s (Q A )yo ~ ln p A / p A for running coupling. 
Thus, when expressed in natural units, the two values are rather close to each other, and 
the fixed coupling estimate could be even slightly larger. 



5.2.2 y > y' c : Proton in the BFKL regime 

So long as 2cby <^ p A (with y > y' c though), the proton saturation momentum remains 
considerably smaller than that of the nucleus, so the proton with k± ~ Q s (A,y) is in the 
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BFKL regime (3.30), where Eq. (3.27) applies. Therefore, the ratio lZ sa t(A, y) is known 
analytically within this whole range, up to an overall constant factor. In order to simplify 
the discussion, we shall study this ratio in particular intervals of rapidity, so that Eq. (3.27) 
reduces to either Eq. (3.36) or Eq. (3.32). 



• The interval p A 2/b < 2cby < p\. 

In this region, tp p is described by Eq. (3.36) with Inz = ln(Q 2 (A,y)/Q 2 s (p,y)) 
T p (y). Since we are in a regime where 2cby ^> p\ we can expand ta(v) = 



ta(v) ~ 



obtain 



\nz = 



p\ 



+ 



2cby + p% to 



(5.19) 



2 V2dy~ 

(Sub leading terms in this expansion, including those arising from the y 1//6 contributions 
to the saturation momentum, cf. Eq. (3.29), result in multiplicative factors in lZ sat (A, y) 
that go to 1 when p A ^> 1.) Then it is straightforward to show that 



lZ sat (A,y) 



A 1 / 3 (2c%) 1 /3 



cxp 



7s 



p\ 



+ 



Pa 



P\ 



2 {2cb y y/ 2 3(2D) 3 / 2 2cby 6(3 (2c%) 3 / 2 



(5.20) 



Notice that with increasing y, the ratio is decreasing. The main difference with respect 
to the corresponding result for fixed coupling, Eq. (5.9), is that for running coupling 
the dominant y-dependencies do not compensate in the ratio Q 2 (A,y)/Q 2 (p,y). Thus, at 
least for not too large y, Eq. (5.20) preserves a relatively fast decrease with y, which is 
dominated by the first term in the exponent, proportional to the "anomalous dimension". 
The other terms in the exponent, due to diffusion effects, are parametrically smaller than 
the first term. However, since these terms are exponentiated, their contribution to TZ sa t 
cannot be neglected in general. 

At the low end of our interval, i.e. when 2cby ~ p X A l ', the ratio becomes 



n s&t {A) 



Pa 



3/10 



A 1 / 3 



exp 



Is 4/5 . 
JPA + 



3/5 



1 



2/5 



3(2£>)3/2^ + 6(3 Pa 



for 2cby 



12/5 

Pa ■ 



(5.21) 



Notice that the very leading behavior A -1 / 3 ~ exp(— p^), has settled in already at this 
rapidity, and what remains in the exponent is fractional powers of p^- These powers 
(4/5, 3/5 and 2/5) are smaller than 1, so that 1Z saX (A) is (much) smaller than 1, but the 
whole exponent is positive so that lZ sa ,t(A) is (much) larger than A^ 1 ^ 3 . This behavior is 
not unexpected though; one would really need to go back to lower rapidities such that 
2cby ~ p A , in order to see the power of A deviate from —1/3 (cf. Eq. (5.20)). 



At the upper end of our interval, i.e. when 2cby ~ p A , the expression becomes considerably 
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simpler since one can keep only the first term in the exponent of Eq. (5.20). One has 



n sat (A) 



exp 



7, 



A 1 ^ 4p~a 

which clearly shows more suppression than Eq. (5.21). 
• The interval p\ < 2cby < p\. 



for 2cby « p\. 



(5.22) 



In this region the proton is in the perturbative geometric scaling region and its gluon 
occupation factor is given by Eq. (3.32). The Cronin ratio becomes 



lZ sat (A,y) 



A 1 / 3 



2 ^/2dby~ 



- -1 


'is p\ 


exp 


2 ^2cby_ 



(5.23) 



It is again clear that lZ sat (A, y) is a decreasing function of y. At the low end of the interval, 
i.e. when 2cby ~ p\, we obtain again Eq. (5.22). On the other hand when we approach 
the upper end p\, we have 



TZ sat (A) oc 



A 1 / 3 



for 2cby ~ p\, 



(5.24) 



where the constant of proportionality in the last equation cannot be determined (since 
the proton approaches a non-linear regime), but it is a number greater than 1. 



5.2.3 2cby > p\ : Proton at saturation 

Even though this region is not analytically under control, since the proton distribution 
becomes non-linear, it is clear that, when further increasing y, the ratio will be (slowly) 
decreasing until it stabilizes, for 2cby ^> p\, at the very small value: 



TZ sat (A, y) ~ for 2cby > p\ , (5.25) 

which is simply the factor introduced by hand in the definition (1.1) of TL v a- Note that 
also the overall normalization is now under control. This reflects the fact that, in this 
limit, the nuclear and proton saturation scales coincide with each other, cf. Eq. (3.25), so 
the corresponding occupation factors at k± ~ Q s (A,y) coincide as well. 

Moreover, by the same argument, it is clear that for 2cby ^> p\, the limiting value 1Z p a = 
l/A 1 / 3 is reached for any k± < Q g (A,y) — this includes the saturation region, and also 
the window (3.30) for the BFKL regime — since in this whole range of momenta the 
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gluon distribution is universal (for sufficiently large y), i.e., it depends upon the specific 
properties of the hadron at hand only via its saturation momentum (cf. Sect. 3.3). This 
will be further discussed in Sect. 6. 

Let us finally compare the asymptotic regimes for fixed and running coupling: By inspec- 
tion of Eqs. (5.11) and (5.25), it becomes clear that the ultimate value of TZ sat is much 
smaller for running coupling, and it is reached after a much longer evolution. Note, how- 
ever, that if one measures the rapidity in units of l/a s , with a s evaluated at the nuclear 
saturation scale Q s (A,y) as natural, then for large y : 



and the asymptotic regime in Eq. (5.25) corresponds to a s y 3> p\. This is parametrically 
of the same order in A as the evolution required to reach the limit (5.11) in the fixed 
coupling case. Still, the presence of the large factor (3 ~ 48 in the fixed coupling condition 
f3a s y 3> p\ makes that, in practice, the lowest value of 7Z sat is reached considerably faster 
with fixed coupling evolution than with the running coupling one. In view of this, the 
running coupling evolution is indeed slower. 



5. 3 The flattening of the Cronin peak 



Since restricted to momenta along the saturation line k± = Q S (A, y), the previous analysis 
could not tell us whether a maximum actually exists or not, neither describe its detailed 
properties, like the position of the peak and its actual shape. To study this, we need the 
nuclear gluon distribution for y > and generic momenta around Q s (A,y), that we shall 
compute in this section to lowest non-trivial order in a s y, that is, after a single step in the 
quantum evolution. Since, as we shall see, the quantum evolution and the twist expansion 
interfere with each other, it is convenient to choose the same "small parameter" , namely 
1/ ' Pa *C 1, to control both of them. If AY = a s y denotes the rapidity increment that 
we shall consider (this satisfies AY <C 1), then our objective in what follows will be to 
construct the nuclear gluon distribution (pA(k±, AY) in such a way that, for AY ~ 1/pa 
and k± < Q S {A, y), all the terms to order 1/pa are correctly included. This means that 
when we use this result to construct the 7£ p A~ratio (and therefore we multiply by pa] see 
e.g. Eq. (2.22)), we just lose contributions that vanish when pa 3> 1. 

This implies, in particular, that only the terms of 0(1) and of 0(1/ pa) need to be kept 
in the initial condition (p in = <pA(k±,Y = 0). (Throughout this subsection, we shall often 
omit the subscript A on the various formulae, to simplify writing.) We thus write: 
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<An = V?o + — <Pi + O (p' A 2 ) , (5.27) 
Pa v 7 

where ipo is the same as the saturating piece in Eq. (2.13) (i.e., the piece involving the 
function T(0, z)), which is of 0(1), while tpi is the term with n = 1 in the twist expansion 
shown in Eq. (2.15) (which, we recall, includes the bremsstrahlung spectrum together 
with higher twist effects to all orders in 1/z) and is of 0(1/ pa)- 

To study the non-linear evolution with y, it is convenient to use the Kovchegov equation 
[33], since this is a closed equation. This can be written both in momentum space, i.e., as 
an equation for <p(k±, Y), and in coordinate space, as an equation for the dipole scattering 
amplitude which appears in Eq. (2.2). In what follows, both forms of this equation will 
be used, for convenience. We start by rewriting Eq. (2.2) as: 

In the following we shall omit the overall factor \/ot s N c , which we can restore at the end. 
In terms of N(r, Y) = J\f xy {Y) (with r± = x± — y± ) , the Kovchegov equation reads: 

dM xy f d 2 z ± 



dY 

where JC xyz is the dipole emission kernel : 



/ 1C xyz {N xz + M zy - M xy - M xz M zy ) , (5.29) 



XVZ ~ (x ± -z ± )*(z ± -y ± y (5 -' U) 

One can transform the above integration over z± into an integration over the sizes of the 
emitted dipoles r\ and r 2 : 

dj\f(r) f dr\ dr 2 2 



r 2 £d£J (e ri )J (£r 2 )J (£r) 

oY J ri r 2 

x [jV(n) + M(r 2 ) - M{r) - JV(ri) N{r 2 )\ . (5.31) 

Through Eq. (5.28), this gives a non-linear equation for tp(k±, Y), in which the structure 
of the non-linear term is actually simpler [33, 38] : 

fp = K®v?-Jy, (5.32) 

where now /C is an operator with the BFKL eigenvalue spectrum, but for our current 
purposes we prefer to view it as a multiple integration determined from Eqs. (5.28) and 
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(5.31), that is : 

K <g> <p = 2 / — —rdred£J (£n) J (£r 2 ) J (£r) J (kr) 



dr\ dr 2 
ri r 2 

x[Af(n)+Af(r 2 )-Af(r)}. (5.33) 



To obtain the change A(p in <p after a small step in rapidity AY <C 1, it is enough to 
iterate the evolution equation only once. This gives: 

Ay? = AY ^/C ® v?i„ - ^fj + O ((AY) 2 ) , (5.34) 

with (fi n the initial condition. To keep the discussion simple, we shall imagine that AY ~ 
1/Pa, which is the amount of rapidity (up to lnp^ factors) that gives a suppression of 
the Cronin peak to a value of order 0(1), due to the fast evolution of the proton (cf. 
Eq. (5.5)). Then we have: 

= AY (iC ® v?o - \ V?o) + O (j^ , (AY) 2 ^ , (5.35) 
and therefore the evolved gluon distribution is finally written as 

<P = V?o + — <Pi + AY^o + o(\, — , (AY) 2 ] , (5.36) 
Pa \p A Pa J 

where we have introduced the shorthand notation 

<Po = /C <g> ipo - - if 2 . (5.37) 

Eq. (5.36) fulfills our objective in that all the neglected terms there become of order 1/ ' p\ 
when AY ~ 1/ Pa (for generic momenta k± ~ Q s (A,y)). 

To compute the action of the BFKL operator on </?o we shall use Eq. (5.33) in which 
M — > Ao = 1 — exp(— r 2 Q 2 /4). Since the three separate integrations (corresponding to 
the three terms in the square bracket) in (5.33) are individually ultraviolet divergent (the 
divergences only cancel in their sum), we shall regulate the r 12 integrations by letting 

dri dr 2 f dr x dr 2 

> hm / - T ^2- -r^r, (5.38) 

r\ r 2 e^o+ J r\ 2e r 2 e 

with the limit to be taken at the end of the calculation. After performing each of the three 
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(multiple) integrations separately, and adding the respective results, we obtain : 

i- 1 r2 ( e ) f z Y 2e 

x 1 - 2 z e T(l - e) $(e, l,-z) + z 2e T(l - 2e) $(2e, 1, -z)] , (5.39) 

where, as usual, we have used the scaled momentum variable z = k 2 /Q 2 {A). Here, 
$(e, 1, — z) =i iq(e, 1, — z) is the confluent hypergeometric function, whose properties 
will be reviewed in the Appendix. The e — > limit in this last expression is finite. After 
also combining this with the non-linear contribution, which simply reads — (1/2)T 2 (0, z), 
we finally find : 

y3o = ^ - \ [ie + In* + r(0, z)f + \ $ 200 (0, 1, -*), (5.40) 

where 7^ = 0.577 is the Euler-Mascheroni constant, and <3> 200 (0, 1, —z) is a particular 
derivative of the confluent hypergeometric function, to be defined in the Appendix. It is 
illuminating to study the small- and large-z limits of <po, which are given by 



¥0 = < 



7T 2 Z 2 

-- T + ... liz «l 



- + ^ H if z > 1. 



(5.41) 



Vz 2z 2 



When z is small, we simply obtain a constant correction. This is consistent with the fact 
that the low-/cj_ limit of tp is not affected by the evolution — this is rather fixed by 
Eq. (5.28) as ip(k, Y) — > In /c 2 when k± — > — , but it is only the scale in this logarithm 
which changes with y : this is the saturation momentum Q S (A, y), since this is the scale 
at which the dipole starts to be completely absorbed: N(r,Y) ~ 1 for r > \/Q s (A,y). 
This discussion suggests that the term (7r 2 /12) AY coming from the one-step evolution 
could be interpreted as the first contribution to the evolution of the saturation scale. 

The crucial point about the correction induced by evolution, as displayed in Eqs. (5.40)- 
(5.41), is that this contains power-law tails at large z, despite the fact that it has been 
generated by evolving the compact (saturating) part of the initial distribution alone. This 
is due to the fact that the (perturbative) gluon interactions — as encoded in the BFKL ker- 
nel — are long-ranged, in both the momentum and the coordinate space (see Eq. (5.30)). 
A similar phenomenon is the generation of power-law tails in impact parameter space 
when solving the Kovchegov equation (5.29), as found in Ref. [52]. 

Thus, in addition to the original "twist" contributions encoded in the function (pi, which 
are parametrically suppressed at large A, the evolved gluon distribution in Eq. (5.36) 
contains other power-law contributions which are not parametrically suppressed, because 
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they are generated by the evolution of the gluons which were originally at saturation: 
ipo and (fa are parametrically of the same order, and the induced piece Atp = AY(p is 
small only as long as AY itself is small. But although we cannot rely on the previous 
approximations to go up to larger values of Y — a s y, it is clear that, for a s y ~ 1 and 
momenta k± ~ Q S (A, y), the twist terms induced by the evolution become as large as the 
saturating distribution. (Incidentally, this justifies the fact that, for a s y ^ 1, we can rely 
on the linear, BFKL, evolution to approach the saturation region from the above; see also 
the discussion after Eq. (3.17)) 

We are now in a position to explain the rapid flattening of the Cronin peak, as seen in the 
numerical simulations in Ref . [25] . We have already noticed, in the analysis of the Cronin 
peak in the MV model in Sect. 2.3, that the presence of long-range tails has the effect to 
flatten the peak, and to push it towards larger values of k±. At y = and for a large value 
of A, a well pronounced peak exists just because the compact contribution to the gluon 
distribution due to the saturated gluons (i.e., ip = (/? sat ) is the dominant contribution at 
k± ^ Qs(A) (see Figs. 2 and Figs. 4). However, with increasing y, the power-law tails are 
enhanced, and the effect of this enhancement is particularly pronounced in the early stages 
of the evolution (a s y <C 1), when the modes at high k± receive large contributions from 
the evolution of the gluons which were originally at saturation. In this way, the exponential 
tail which was present in the saturating distribution at y = is rapidly washed out, and 
replaced by a power law tail which starts already at k± ~ Q s (A,y). When this happens, 
the peak in the ratio TZ p a flattens out completely and disappears. In view of the previous 
estimates, we expect a complete flattening for a s y ~ 1. But the effect shows up already 
for smaller rapidities, where our one-step evolution is a good approximation. 

To illustrate this, we have represented in Fig. 7 the ratio 1Z p a (for pa = 6) obtained by 
using the result in Eq. (5.36) for the nuclear gluon distribution together with the DLA 
approximation for the proton distribution 14 . For comparison, we have also represented on 
the same plot the ratio which is obtained when the non-evolved, MV model, distribution 
is used for the nucleus. The rapid suppression of the peak, due to the fast rise in the proton 
distribution, is clearly seen in both cases. But in the absence of nuclear evolution the peak 
is always there; just its amplitude gets smaller and smaller. By contrast, when using the 
evolved nucleus distribution from Eq. (5.36) the flattening of the peak is manifest, and in 
fact the maximum has almost disappeared after an evolution AY = 2/p A ~ 0.3. 

It is also visible on Fig. 7 that with increasing y, the position of the peak moves towards 
larger momenta — in agreement with our expectation that the peak should follow the 



To be able to perform this calculation also for very small values of Y, where the saddle point 
approximation is not appropriate anymore, we have used the exact solution to the DLA equation, 

which is well known to read (compare to Eq. (3.15)) : f p (k±,y) = -^-jj- pf ^o(\/4a s y/o), with 
p = ln(kj_/Qp) and Iq(z) the modified Bessel function of the first kind. 



66 




Fig. 7. The Cronin ratio 1Z p a{z) below and near the saturation scale for pa = 6. The black 
(solid) lines correspond to an evolved nuclear wavefunction by Ay <C 1. The red (dotted) lines 
correspond to an unevolved one (MV) . The proton wavefunction is always given by the full DL A 
solution. The curves, from top to bottom, correspond to Ay = n/(2pj\) with n = 0, 1, ...,4. 

nuclear saturation momentum — , but it does that only very slowly. This is so slow because, 
first, the nucleus evolves only little for such a small rapidity increment, and, second, there 
is an opposite effect due to the DLA evolution of the proton, which for the small values of Y 
of interest here is almost compensating the evolution towards larger k± due to the nucleus. 
Using Eq. (5.36) for the nuclear distribution, together with the DLA approximation (the 
full Bessel function) for the proton, we have been able to compute analytically the position 
of the peak for very small values of y. The result reads (compare to Eq. (2.23)) : 

zq(Y) « 0.435 + + 0.862F - 0.769a —, (5.42) 

Pa V Pa 

where the neglected terms are such that, when Y ~ I/pa, they are all of 0(p~^ 2 ). 

To summarize the previous considerations, the ratio (1.1) can be suggestively rewritten 
as follows: 



$A(k ± ,y) 



(fA{k±,y) 
<PA(k±,y = 0) 



%{k±,y) 



fcpA(k ± ,y = 0) 

= <P P (h±,y) 
<P P (k±,y = 0) 



(5.43) 
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The ratios $a and $ p , which describe the relative evolutions of the nucleus and of the 
proton, respectively, play very different roles for the evolution at small y and momenta 
around Q s (A,y) : $ p rises very fast with y, and is responsible for the rapid suppression 
of TZpA at generic momenta. However, this factor varies only slowly with k±, and thus 
cannot be responsible for the flattening of the peak. $^4, on the other hand, rises only 
slowly with y, but its evolution is quite asymmetric around Q S (A, y) : the (relative) 
evolution is substantially larger at momenta above Q s (A,y) than below it, and because 
of this asymmetry, the peak gets progressively tilted (its 'upper' side at larger k± rises 
faster than the lower one) and also flattens, until it eventually disappears. 



6 High-A;^ suppression in the nuclear gluon distribution at small x 

In this section we shall extend our analysis of the evolution of the ratio 7Z p A{k±,y) with 
increasing y to arbitrary values of the transverse momentum k±. The various regimes of 
evolution in the y — In k\ plane are illustrated in Figs. 1 and 10 for fixed and running cou- 
pling, respectively. In Sect. 5, we have followed the evolution along the nuclear saturation 
line k± = Q s (A,y). Below, we shall consider all the physical regimes exhibited in Figs. 
1 and 10, but our main focus will be on the high momentum region at kj_ ^> Q s (A,y), 
which is the most interesting region for the phenomenon of "high-A;_i_ suppression" . The 
natural way to follow this evolution is to increase simultaneously y and k±, in such a way 
that the ratio z = kj_/Q^.(A,y) remains constant. This ensures that, with increasing y, 
the nucleus remains in the same physical regime, which could be either the saturation 
regime (for z < 1), or the (nuclear) BFKL regime (for z > 1). 

In performing this analysis, it is useful to keep in mind the limitations of the analytic 
approximations that we shall use for the gluon occupation factor: The transition regions 
from saturation to BFKL, and also from BFKL to DLA, are not analytically under con- 
trol, nor are the very early stages of the evolution. For the BFKL regime, as defined by 
Eq. (3.20), the saddle point approximation is justified provided d s y > 1, while for the 
DLA regime at p ^> a s y, the corresponding condition reads a s yp ^> 1. We recall that 
p = p(A,k±) = In k\/Q 2 s (A) is much larger for the proton than for the nucleus (for the 
same value of k±). Because of that, we shall be able to study the early (a s y <C 1) evolution 
of the proton (as described by DLA), but not also that of the nucleus (unless k± is so 
large that the nucleus itself enters the DLA regime). 

The analysis below will not only confirm the general arguments in Sect. 4 about the 
variation of the ratio 7Z p a with y, k± and A, but will also allow us to explicitly compute 
this ratio and, in particular, its limiting values at large y and/or large k±. In the fixed 
coupling case, and for the linear regime, our results below will confirm and extend some 
of the conclusions originally obtained in Refs. [16,24], with which the present analysis 
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has some overlap. Also, still for fixed coupling, our results will be seen to be consistent 
(within their limited range of validity and to the present accuracy) with the numerical 
analysis in Ref. [25] based on the Kovchegov equation. 

6.1 Fixed coupling 

In this section we shall explore the various regimes for the fixed coupling evolution exhib- 
ited in Fig. 1, and study the variation of the ratio 7Z p A(k±,y) with y and k±, as well as 
its parametric dependence upon A. 

6.1.1 k± < Q s (A,y) : The nuclear saturation region 

The nuclear distribution at saturation is not explicitly known (except for very small 
momenta k± <C Q S (A, y), where Eq. (3.1) applies), but this is some slowly varying function 
of z = k\/Q 2 s {A,y) (because of geometric scaling), which moreover is simply a constant 
when the evolution is performed along a line parallel to the saturation line. As for the 
proton, it can be either at saturation, or in the linear regime above saturation, depending 
upon the value of k±. 

(I) k± < Q s (p,y) ■ Proton at saturation 

In this regime, the gluon distributions in both the proton and the nucleus are universal 
functions of k\/Q 2 s {A, y) (with Q s (A,y) — > Q s (p,y) in the case of the proton), which are 
only slowly varying. Thus, clearly, the ratio (1.1) is of order 1/A 1 ^ 3 , with only weak, loga- 
rithmic, dependencies upon A and k±. This becomes explicit for sufficiently low momenta, 
k± "C Q s {Pi y)) where one can use Eq. (3.1) for both the proton and the nucleus, and thus 
find that TZ p A{k±,y) is monotonously increasing with k±, and slowly, but monotonously, 
decreasing with y at fixed k±. In particular, along the proton saturation line k± = Q s (p, y), 
which defines the upper boundary of this region, one finds: 

K P A(k_L = Q s (p,y),y) ~ (independent of y) . (6.1) 

(II) Q s (p,y) < k± < Q s (A,y) : Proton in the linear regime 

The corresponding analysis is very similar to that of the evolution along the nuclear 
saturation line, as discussed in Sect. 5. For y < y c , cf. Eq. (5.1), the proton can be either 
in the DLA, or in the BFKL, regime (depending upon the value of k±), while for larger 
y, it is always in the BFKL regime (see Fig. 1). 
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A straightforward analysis shows that the ratio is rapidly decreasing with y at any fixed 
k±, due to the fast evolution of the proton. It is also interesting to consider the large 
y behavior at fixed z. Then, the proton is in the BFKL regime, so using Eq. (3.17) 
(with Q S (A, y) — * Qs(p,y)) and writing <PA(k±,y) = (l/a s N c )<p (z), one obtains (with 
7 = 7 S ~ 0.63 from now on) : 



>o a \ Pa z J( Po(z) j (p A + In z) 

K pA {k ± ,y) ~ -rrj exp ^ -— — 

A— Pa - In l/z + A { 2(3a s y 

. Pa z^ (z) 2 



Pa - In ' 



for 2(3a s y^>p 2 A , (6.2) 



(in this regime pa > hal/z), which is indeed consistent with the previous result (5.11) on 
the nuclear saturation line z — 1. 



6.1.2 Q s (A,y) < k± < Q g (A,y) : The nuclear BFKL region 

For y < y c , the proton is in the DLA regime for any k±, while for y > y c , it can be either 
in the BFKL, or in the DLA regime, according to the value of k±. 

(I) y < y c : Proton in the DLA regime 

In this region, the ratio TZ p a is formed by using Eq. (3.17) for the nucleus and Eq. (3.15) 
for the proton. After straightforward manipulations, similar to those leading to Eq. (5.2), 
the result can be cast into the form: 

TZ pA (z,Y) ~ Paz 1 - 1 (ln^ + A) exp jcF- ^(ln^ + cF + p A ) - ^| , (6.3) 

where z = k\/Q 2 s {A, y) and Y = a s y. Note that the overall normalization is not under 
control. Eq. (6.3) applies for k± in the nuclear scaling window (3.20), that is: 



< In z < In 



Q 2 (A, y)/Q 2 M, V) = "cY, and 1 < Y < Y c = p A /vc. (6.4) 



It can be checked that, within this physical range 15 the function TZ p a(z, Y) is monotonously 
increasing with z, and also rapidly decreasing with Y at any fixed z. This behavior is in 
agreement with the general results in Sect. 4, and is illustrated in Fig. 8. 

(II) y > y c & Q s (A,y) < k± < Q g (p,y) : Double-scaling regime 



15 The function in Eq. (6.3) develops a maximum at very large z, well outside the window (6.4). 
This is, of course, unphysical. 
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Fig. 8. The ratio 1Z p a(z,Y) in Eq. (6.3) as a function of z = k\/Q 2 s (A,y) for pa = 6 (i.e., 
Y c f» 0.75) at different rapidities. The lines are for Y = 0.1 + 0.15 n, with n = 0, 1, ... ,4, from 
top to bottom. The bottom line corresponds to Y = 0.7 ~ Y c . The ratio 7£ p a is measured in 
arbitrary units, since the normalization of Eq. (6.3) is not under control. Note the rapid decrease 
with increasing Y in the early stages of the evolution. 

This region is interesting in that the proton and the nucleus are both in the same physical 
regime, namely in the window (3.20) for (approximate) geometric scaling, so they are 
described by the same approximation — the BFKL expression (3.17) — , and the nor- 
malization ambiguities cancel in the ratio 7Z p a- Thus, in this regime, one can also predict 
the amplitude of this ratio, and not only its functional dependencies. Besides, this is also 
the regime for which the quantum evolution above the saturation line has been originally 
invoked, by Kharzeev, Levin and McLerran [16], as a possible mechanism to explain the 
high-fcx suppression observed in the RHIC data. 

Using Eq. (3.17) for both the nucleus and the proton, one finds : 



pi 



\nz + A 



A-r 1 lnz + p A + A 



cxp 



Pa 
2(3Y 



;21n z + p A ) 



(6.5) 



where we have used \n[kj_/Q 2 (p,y)} ~ \nz + Pa- Eq. (6.5) is valid within the range: 

< \nz < vcY — pa and Y > Y c = pa/vc. (6.6) 



It is easily checked that the ratio (6.5) is an increasing function of z for arbitrary Y, and 
a decreasing function of Y for arbitrary z > 1. This behavior is illustrated in Fig. 9. 
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A noticeable difference with respect to the previous case, cf. Eq. (6.3), is that, in this 
double-scaling regime, the dominant y-dependencies of the nuclear and proton distri- 
butions, which are exponential, have cancelled in the ratio 1Z p a- Indeed, according to 
Eq. (3.17), these dependencies are encoded in the respective saturation momenta, which 
however evolve in the same with y, cf. Eq. (3.14), and thus their ratio Q 2 S (A, y)/Q 2 s (p, y) = 
Q 2 (A)/Qp ~ A x l z p A is independent of y. 

The residual dependence on y in Eq. (6.5) comes from the diffusion term in Eq. (3.17), 
and is rather weak. (This is also in agreement with the discussion in Sect. 4; see especially 
Eq. (4.10).) Thus, for Y > Y c , the suppression slows down, and for sufficiently large Y 
(such that the exponent in the second line of Eq. (6.5) become negligible), 7Z p a stabilizes 
at a small value which depends slowly on z: 

It can be verified that most of this suppression has been achieved already in the early 
stages of the evolution, while the proton was still in the DLA regime. 

When approaching the nuclear saturation line from the above, Eq. (6.5) yields: 

n - iz lX) * (Psp e * ^ (6 ' 8) 

As expected, this is the same function of A and y as obtained when approaching the 
saturation line from the below, cf. Eq. (5.9), with the noticeable difference that, now, also 
the normalization of this result is under control. 

It is finally interesting to evaluate the ratio (6.5) on the proton geometric line (i.e., k± = 
Qg(j), y), or \nz = ucY — pa), which is the upper boundary for the double-scaling regime 
considered here. In the large Y limit, one obtains: 

\\mK pA (k ± ~Q g ( p ,y),y) = { J 3 A pA)S , 5 = 1-7-^0.21. (6.9) 

The dominant power of 1/A 1 ^ 3 is lower on this upper boundary than on the lower one 
(5 ~ 0.21 as compared to 1 — 7 ~ 0.37), since k± is increasing faster along the geometric 
scaling line, and 7Z p a is increasing with k±. 

(Ill) y > y c h Q g (p,y) < k± < Q g {A,y) : Proton in the DLA regime 

In this regime, the ratio 1Z p a(z, Y) is given again by Eq. (6.3), but which now applies to 
a different kinematical range : ucY — pa < In 2; < ucY and Y > Y c . 
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Fig. 9. The ratio 1Z p a(z,Y) in the double scaling window for p A = 6 as a function of z at 
different Y > Y c « 0.75. The solid lines correspond to Y = 0.75 + 0.3 n, with n = 0, . . . ,4 from 
top to bottom. The dashed line is the asymptotic (Y — ► oo) profile in Eq. (6.7). The upper limit 
of the double scaling window z max = exp{z^c(Y — Y c )} rapidly increases with increasing Y. For 



instance, 



11 for Y = 1.05 (n = 1), and z n 



127 for Y = 1.35 (n = 2). 



i.<? A;_|_ > Q 9 (A, y) : The nuclear DLA region 

This is the high momentum limit, in which both the proton and nucleus are in the DLA (or, 
more properly, DGLAP) regime, so the ratio can be again computed without uncertainties 
related to the normalization. One finds: 



TZ pA (z,Y) = exp 



4Y(\nz + cY) - J4Y(\nz + cY + p A ) 



(6.10) 



valid for \nz > vcY . It is clear that the extra term p A in the exponential comes from the 
difference between two saturation scales. It is easy to check analytically that the ratio is 
an increasing function of z for any Y, but a decreasing function of Y for any fixed z. In 
particular, at very large z such that In 2; ^> cY + p A , one obtains 



K pA ~ e 



-PA 



(6.11) 



Therefore, when z — > 00 the ratio approaches one from below. 
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In k 



Fig. 10. Various regimes of evolution for running coupling. The boundaries which separate the 
three physical regimes (CGC, BFKL, and DLA) are shown by the red (thin) and blue (thick) 
curves for the proton and the nucleus, respectively. For each color (thickness), the left curve 
corresponds to the saturation line (k± = Q s (p, y) or k± = Q S (A, y)), while the right curve is the 
boundary between BFKL and DLA regimes (which is less under control). The shaded bands on 
the right of the saturation lines are the respective geometric scaling regimes, which overlap with 
each other in the "double scaling" regime (the area with a darker shade). 

6.2 Running coupling 

The various physical regimes for evolution with running coupling are exhibited in Fig. 10, 
which should be compared to the corresponding map for fixed coupling, in Fig. 1. The 
differences between the two cases come from several sources: 

First, the two saturation scales evolve differently with y (that for the proton growing 
faster) and eventually converge towards a common limit at very large y (cf. Eqs. (3.23)- 
(3.25)). A similar behavior holds for the corresponding lines for geometric scaling, and also 
for the boundaries separating BFKL from DLA behavior. This implies that, for sufficiently 
large y, the gluon distributions in the proton and in the nucleus become identical for all 
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the interesting momenta. (Differences between the two distributions persist only in the 
(common) DLA regime, which with increasing y is rapidly pushed towards asymptotic 
momenta.) This leads to the important conclusion that the ratio (1.1) approaches the 
universal constant 1/A 1 ^ 3 at very large y, for all non-asymptotic momenta. 

Second, in the case of a running coupling, the window for geometric scaling above satura- 
tion becomes narrower (cf. Eq. (3.33)), and, in contrast to the fixed coupling case, it does 
not coincide with the BFKL regime anymore : the latter is now defined by Eq. (3.30), 
and includes a domain where scaling violations are important. In Fig. 10, the regions 
of extended geometric scaling are the shaded bands lying on the right hand side of the 
nuclear and proton saturation lines, respectively. 

Third, the precise boundaries of the BFKL and DLA regimes, and also the transition 
between these two regimes, are not under control. Whereas we expect the BFKL approx- 
imation (3.27) to be valid for < Inz < Cir^ /3 (y) (with z = k 2 ± /Q 2 s (A,y)), and the DLA 
approximation (3.26) should apply at In z > C2TA.(y), the constant factors 0^2 in these 
conditions are not determined, and we do not dispose of analytic approximations for the 
intermediate region at ciT^ 3 (y) < \nz < C2Ta(v)- 

To avoid cumbersome discussions and ambiguities associated with the location of the bor- 
derlines for the various regimes, in what follows we shall focus mainly on the situations 
in which both the proton and the nucleus find themselves in the same physical regime 
(which could be CGC, BFKL, or DLA). These situations are the most interesting ones, 
as they allow for a calculation of the ratio 1Z p a which is free of renormalization ambigu- 
ities. Besides, they are the only ones to survive at large y: The other regimes visible on 
Fig. 10 are squeezed in between the proton and nucleus saturation lines, or in between 
the respective borderlines for BFKL behavior, and therefore shrinks to zero in the high 
energy limit. 

The analysis to follow is particularly important because, with running coupling, we have 
no general arguments to understand the behavior of the ratio 7Z p a, as we had for a fixed 
coupling. As we will see, in spite of interesting differences, the global behavior of 1Z p a as 
a function of y, k± and A is the same as in the fixed coupling case discussed in Sect. 4. 

6.2.1 The double CGC region 

On the left of the proton saturation line, both the proton and the nucleus are in the color 
glass condensate regime. For sufficiently low momenta k± <C Q s (p, y), one can compute the 
ratio TZpA by using Eq. (3.2) for both hadrons. Then one can easily check that TZ p A{k±, y) 
has the same behavior as in the corresponding regime at fixed coupling. It is further 
interesting to follow the proton saturation line k± = Q s (p,y), where Eq. (3.4) holds for 
the proton. One finds (see Eq. (3.23) for notations) 
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where we have added a constant k to the nuclear distribution in Eq. (3.2) in such a way to 
match the boundary condition (3.4) when k± — > Q s (A,y). Of course, this is only a crude 
approximation — the approach towards the saturation line is not under control — , but 
this is interesting as it allows us to qualitatively study the behavior at large y, where the 
two saturation momenta converge to each other. For 2cby ^> p\, one finds: 



K^ = QMy),y)~^(i + a ^^), (6.13) 

which decreases with y, and for 2cby ^> p\ approaches the universal value I /A 1 / 3 , as 
expected. 



6.2.2 The nuclear BFKL region 

Here we consider the region in between the two blue (thick) curves in Fig. 10. For relatively 
small rapidities y < y c with 2cby c ~ p\, the proton is in the DLA regime for any k± 
within this region. (See the discussion at the beginning of Sect. 5.2 for more details on the 
separating rapidities.) For larger rapidities, the proton BFKL regime starts to appear, but 
this is theoretically under control only for y > y' c , with 2cby' c ~ p\ 2 ^ . The most interesting 
"double scaling" regime appears at the even larger rapidity 2cby" c ~ p\. 

(I) y < y c : Proton in the DLA regime 

Straightforward manipulations using Eq. (3.27) for the nucleus and Eq. (3.26) for the 
proton yield (with z = k\j 'Q 2 s (A,y), as usual) 



n pA (z, y)~(^lnz + r A (y)^j ^(y) Ai U + 



lnz + A_\ ^ 



DrTiy) 



x exp |r A (y) - p A - \/^by \n(\nz + r A (y)) - | , (6-14) 

where the normalization is not under control. This is valid provided the nucleus is in the 
BFKL window (3.30), which implies (ci is a constant of 0(1)) : 



< \nz < ^T^iy). (6.15) 
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It can be checked that, within this physical range, the ratio 1Z p A(z,y) is monotonously 
increasing 16 with z, and also rapidly decreasing with y at fixed z. 

In particular, when approaching the nuclear saturation scale from the above, one recovers 
the same parametric dependencies upon A and y as in Eq. (5.16). 

(II) y > y' c : The double scaling regime 

As before, this regime is interesting since the uncontrollable normalization factors cancel 
out in the ratio between the gluon distributions for the nucleus and the proton. Using 
Eq. (3.27), one finds : 



lZp A (z,y) ~ 



ta 



1/3 



Ai ( {l + ^#) 



A 1 / 3 



Ai6 + 



In z+ta~t p + A 



Dt, 



1/3 



x exp 



i(t a ~ T p ) + 



2(\nz + t a -t p ) 2 2\n 2 z 



(6.16) 



which is free of normalization ambiguities. Since now y is relatively large, one can replace 
in this and the subsequent equations Ta — t p ps p 2 A /2^2cby and r p ps \j2cby. 

The upper limit for the validity range of this expression is fixed by the borderline for the 
proton BFKL regime (see Fig. 10). Therefore, Eq. (6.16) is valid for: 



< lnz < ci r 2 J\y) - (r A (y) - r p (y)). 



(6.17) 



It can be checked that, the ratio (6.16) is an increasing function of z, and a decreasing 
function of Y. This can be explicitly seen on the numerical plot of the ratio in Fig. 11. 
Note the rapid decrease of the ratio with increasing y. As compared to the corresponding 
plot for fixed coupling, in Fig. 9, one can see that with a running coupling the decrease 
of TZ p a(z, y) is pursued up to larger values of y, until the ratio eventually stabilizes at the 
^-independent (and comparatively smaller) value 1/A 1 ^ 3 . 

This evolution can be more explicitly studied by using the approximate form of Eq. (6.16) 
valid in the "double scaling" window, which is the regime where both the nucleus and 
the proton exhibit geometric scaling. This window opens at a rapidity 2cby" c ~ p% and 
extends over the following interval in z (c^ is a constant of 0(1)); see Eq. (3.33)): 



< In z < In 



Q 2 Jp,y)/Q 2 s (A,y)] = csr^iy) - (TA(y)-r p (y)), 



(6.18) 



16 Although the ratio (6.14) develops a maximum at large z, it is well outside the physical 
region, similarly to the fixed coupling case. In fact, one can roughly estimate the position of the 
maximum to be lnz ~ aTA(y) with a ~ 0(1). 
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Fig. 11. The ratio lZ pA (z, y) in the double BFKL regime for p A = 6 as a function of z at different 
values of y. The lines correspond to y = 1.1 + 0.3 n, with n = 0, . . . ,4 from top to bottom. The 
third line (n = 2) corresponds to y = y" = 1.7. The range in z shown in the figure is in the 
physical range for all the lines, except for the top one. Indeed, the upper boundary of the physical 
regime is z max = 22.8 for y = 1.1, z max = 50 for y = 1.4, z max = 97 for y = y" c , etc. 

which is represented by the dark shaded area in Fig. 10. In this range, the ratio takes the 
simpler form (cf. Eq. (3.32)) : 



A 



A 1 / 3 \nz + p 2 A /2^/2cT>y + A 



(6.19) 



which when approaching the saturation line from above (lnz — > 0) is consistent with the 
previous estimate in Eq. (5.23). 

For large y, such that 2cby ^> p\, the ratio (6.19) becomes independent of z, and ap- 
proaches the universal limit 1/A 1 ^ 3 , as anticipated. This is also the limit of the full ex- 
pression (6.16) for 2cby ^> p\ and any z in the range (6.17). 



6.2.3 The nuclear DLA region 

By using Eq. (3.26) for both the proton and the nucleus, one easily obtains: 



Tl pA (z,y)=exp ^4&y{ln(ln2 + r A (y)) - \np A } - yj 4by {ln(ln z + r A (y)) - \np p } 



(6.20) 
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which is valid, roughly, for Inz ^ ta{u)- It is easy to check that the ratio is an increasing 
function of z for any y. We have numerically verified that, for In 2; ^ ta(v) at least, the 
ratio is an decreasing function of y. 

At very large z such that In 2; 7a (y), one obtains 

TV ~ e - ln ^^ 62//lnln2 . (6.21) 

Thus, when z — > 00 the ratio approaches one from the below (but only very slowly). This 
concludes our analysis. 
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A Appendix 



Here we study in more detail the gluon occupation factor ipA and the integrated gluon 
distribution function Qa in the McLerran- Venugopalan model, both for fixed and running 
coupling. We first deal with the running coupling case, which turns out to be simpler, 
and then we continue with the fixed coupling one. At the end of the appendix we present 
some useful properties of the confluent hypergeo metric function. Note that, in order to 
simplify writing, we shall often use different normalizations for the various quantities as 
compared to the main text. The differences will be indicated at the appropriate places. 
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The Running Coupling Case 



The twist contribution to the gluon occupation factor, as given in Eq. (2.17), and with 
the constant factor l/(b N c ) omitted since it can be restored at the end of the calculation, 
reads 

^\ z ) =JdtJ (Vtet) In - . (A.l) 



In principle, one should integrate up to ~ t max = e PA , which corresponds to a dipole size ~ 
r mffi = 2/A. However, it is harmless to extend the integration to infinity, so long as we are 
interested in momenta k 3> A z 3> e~ PA . We can do the above integration by replacing 
ln(l/t) — > t~ a . Then we obtain an expression involving the confluent hypergeometric 
function iiq(— a, 1, —z) = <&( — a, 1, —z). Finally we take the derivative with respect to a 
at a = to recover the quantity in Eq. (A.l), and we obtain 

^twist = _ 1 ln 2 z _ 7 ln z + 7 r(0? ^ + g _ f + 1 $ 200 (0) x _^ (A 2) 

where 7 = 7^ = 0.577 denotes the Euler-Mascheroni constant throughout this Appendix. 
In the small-z limit the hypergeometric term vanishes and one is left with 

^twist = _ 1 ln 2 z _ 27lnz + g _ 3f + ^ (A 3) 

Notice that this is negative, but it never becomes larger in magnitude than the contribu- 
tion of the saturation term ip 8 ^ = PaT(0, z). Even at k = A z = e~ PA , one has 

<pT* * ~\p\ ^ AfT- (A.4) 
and the gluon occupation factor is positive and well defined. 

In the large- ,2 region, we can use Eq. (A. 21), given at the end of the appendix, to derive 
the asymptotic expansion of the index-differentiated hypergeometric function. This will 
have the form given in Eq.(A.23). It is straightforward to show that 

$ 200 (0, 1, -z) = ln 2 z + 2 7 ln z + 7 2 - — + 2 V + 0( e - z ). (A.5) 

6 n=1 n z n 

Then, Eqs. (A. 2) and (A.5) imply that all the logarithmically divergent and constant 
terms in <pa cancel, and as expected only the inverse powers of z survive. Therefore one 
has 

„ twist v r(n) 1 Q 2 S (A) I Q^A) 2Ql{A) 
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which exhibits the correct bremsstrahlung spectrum, accompanied by power corrections 
which are all positive. 



Now we turn our attention to the integrated gluon distribution function, which we con- 
veniently normalize as 

«m^- (A - 7) 

For the saturation term of the gluon occupation factor the integration is trivial, and one 
obtains 

QT =[l-e' z + ZT (0, Z)] p A . (A.8) 

In the above we have defined the scaled momentum variable Z = Q 2 /Q 2 (A). For the 
twist contribution as given in (A.l), one can perform first the integration over k, which 
simply replaces J (y/4zt) by \J~Z J\ {yj Azt^j j\ft. Then the integration over t can be done 
as before and gives 



/i twist 
y A 



7T 2 7 2 



- Z In 2 Z + (1 - 7) Z In Z + — + 7 - 1 \Z 



12 



(A.9) 



1 



+ 7 [l - e~ z + Z T (0, Z)\ + - Z $ 200 (0 , 2, -Z). 
When Q 2 » Q 2 (A) «Z>1, the last two equations lead to 

g A = PA + \nZ + 2 1 -j:^± + 0(e- z ), 

n=1 n + 1 z 



(A.10) 



where the first term comes from Q\ &t and the remaining from Q A mst . One can rewrite the 
above equation, in a more illuminating way, as 



Q 2 



1Q 2 (A) IQtiA) 



gA = ln V +2 ^-2 Q 2 3 Q* 



+ 



(A.ll) 



which shows that the higher twist corrections give a small negative contribution (shad- 
owing) to the integrated gluon distribution function. 



The Fixed Coupling Case 



The twist contribution to the gluon occupation factor as given in Eq. (2.15) and multiplied 
by the constant factor p A a s N c , reads 



".max j.n—1 l n ri + 

dtJ (VAzt) e~ l —y^r- (A.12) 
n =i n - Pa 
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The integration in the above equation needs to be cut at ~ t max = e PA which corresponds 
to a dipole size ~ r max = 2/A. However, and in the "worst case" that z is small, each 
term gets its most contribution from the region t ~ n, with the "magnitude" of the n-th 
term being ~ (\nn/p A ) n . Therefore the integration in the terms with n < n max = e pA can 
be extended to infinity, while terms with n > n max do not contribute, since their peak lies 
outside the region of integration. Thus, the cutoff in the MV model translates to a cutoff 
in the series in Eq. (A. 12). Furthermore, this (asymptotic) series is rapidly convergent 
for pa 1, since lnn/p^ < 1 for n < e PA . Thus, it is enough to keep a few terms in 
(A. 12). They can be calculated by the same method we used in the running coupling case. 
The n-th term, apart form the overall suppressing factor p^ n+1 , becomes constant when 
z <C 1, and falls as z~ n (plus all possible higher inverse powers) for z ^> 1. The expressions 
become more and more complicated as n increases, and therefore we shall present only 
the first two terms. We have 



¥a = 

Pa 



* 100 (i,i,-*), 



1 



f(l- 7 )$ 100 (2,l,-^) + ^ 200 (2,l,-z) 



(A.13) 



(A. 14) 



Since $ (a, 1, — z) vanishes at small-z, the limiting value of ip 1 ^ 1 is finite. One easily 
finds that 




(A.15) 



In the large-z limit, we need to expand the hypergeometric functions and we obtain 



tn ( twist Q 2 S (A) Q 4 S (A) 



Q 2 



Q 4 



1 + — f ln tJt-TT + 2 7-2 

PA V Q 2 9 (A) 



+ 



(A.16) 



where both the n — 1 and n = 2 terms contribute to the 1/z 2 = Q^(A)/Q 4 term. This 
correction is positive and the same will be true for all the l/z n ones. 



Now let us calculate the integrated gluon distribution function. The saturation contribu- 
tion is the same as the one we found in the running coupling case and is given in Eq. (A. 8). 
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For the twist contributions we shall again present only the first two terms, which are 



r>n=2 

yA = 



Pa 



z 



Z$ 100 (l,2,-Z) 

2 „,2 \ 



(A.17) 



7T 

12 + ~2 



7 e 



1 



+ (1 - 7 ) $ 100 (2, 2, -Z) + - $ 200 (2, 2, -Z) 

2 



(A.18) 



and similar, but more complicated expressions can also be obtained for n > 3. Each term 
vanishes as Z becomes small, and all terms with n > 2 approach as Z becomes large. 
Keeping the n = 1,2 terms, we can obtain the first power law corrections in Qa in the 
high momentum limit Q 2 S {A) ^> Q 2 . We obtain 



^ = ln^ + 2 7 -^ 



i+ ^( ,n & +a7 - 1 ! 



(A.19) 



As in the running coupling case the higher twist terms give a small negative contribution 
to the integrated gluon distribution function. 



The Confluent Hypergeometric Function <3>(a, b, —z) 



Here we give some useful properties of the confluent hypergeometric function i-F\(a, b, —z) = 
$(a, 6, — z). More precisely we examine its n-th derivative with respect to the first argu- 
ment, that is $ n00 (a,6, —z). The range of parameters that we are interested in is a > 0, 
b > 1 and z > 0. 



The confluent hypergeometric function is defined as 



a 1 a(a + 1) 2 

K^-z) = i-j, + ilF5 ^- 

and its large- z asymptotic expansion reads 



(A.20) 



x 



r(6 - a) 

a(a- 6+1) 1 a(a + l)(a-b + l)(a-b + 2) 
+ z + 2! +" 



(A.21) 



where terms that fall as e 2 have been neglected. 
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We can find the behavior of <3> ri00 (a, 6, —z) in the small- z region from Eq. (A. 20) and in 
the large- z region from Eq. (A. 21). We obtain 



$ n00 (a,b,-z) 



r(6) 



T(b + n) 



z n + 0{z 



n+l\ 



for z <C 1, 



(A.22) 



which vanishes when n > 1, and 

$ n00 (a,6, -z) = z- a 



E 

fe>0 



Ck 



E ln * 



.fe=0 



for z >> 1, 



(A.23) 



where the coefficients and c4 depend on the set of parameters (n, a, b) and can be 
determined from Eq. (A. 21). The function $ n00 (a,6, —z) diverges for large- z as a power 
of In z when a = 0, while it approaches zero when a > 0. 



References 



[1] "Quark Matter 2002", Proceedings of the 16th International Conference on Ultra-Relativistic 
Nucleus-Nucleus Collisions, (Nantes, France 18-24 July, 2002), Nucl. Phys. A 715 (2003). 

[2] See for example, D. d'Enterria, "Hard scattering at RHIC: Experimental review, " arXiv:nucl- 
ex/0309015 and references therein. 

[3] B. B. Back et al. [PHOBOS Collaboration], Phys. Rev. Lett. 91 (2003) 072302 [arXiv:nucl- 
ex/0306025]; S. S. Adler et al. [PHENIX Collaboration], Phys. Rev. Lett. 91 (2003) 072303 
[arXiv:nucl-ex/0306021]; J. Adams et al. [STAR Collaboration], Phys. Rev. Lett. 91 (2003) 
072304 [arXiv:nucl-ex/0306024]; I. Arsene et al. [BRAHMS Collaboration], Phys. Rev. Lett. 
91 (2003) 072305 [arXiv:nucl-ex/0307003]. 

[4] I. Arsene [BRAHMS Collaboration], "On the evolution of the nuclear modification factors with 
rapidity and centrality in d+Au collisions at yJs^N = 200 GeV," arXiv:nucl-ex/0403005. 

[5] R. Debbe and M. Murray, for the BRAHMS Collaboration, talks given at QM2004, Oakland, 
USA, January 2004. 

[6] T. Frawley and M. X. Liu, [PHENIX Collaboration], P. Steinberg and R. Nouicer, [PHOBOS 
Collaboration], and L. Barnby [STAR Collaboration], talks given at QM2004, Oakland, USA, 
January 2004. 

[7] L.V. Gribov, E.M. Levin, and M.G. Ryskin, Phys. Rept. 100 (1983) 1. 

[8] A.H. Mueller and J. Qiu, Nucl. Phys. B268 (1986) 427. 

[9] J.-P. Blaizot and A. H. Mueller, Nucl. Phys. B289 (1987) 847. 



84 



[10] L. McLerran and R. Venugopalan, Phys. Rev. D49 (1994) 2233; ibid. 49 (1994) 3352; ibid. 
50 (1994) 2225. 

[11] A. H. Mueller, Parton Saturation-An Overview, hep-ph/01 11244. Published in QCD 
Perspectives on Hot and Dense Matter, Eds. J.-P. Blaizot and E. Iancu, NATO Science 
Series, Kluwer, 2002. 

[12] E. Iancu, A. Leonidov and L. McLerran, The Colour Glass Condensate: An Introduction, 
hep-ph/0202270. Published in QCD Perspectives on Hot and Dense Matter, Eds. J.-P. Blaizot 
and E. Iancu, NATO Science Series, Kluwer, 2002; 

E. Iancu and R. Venugopalan, The Color Glass Condensate and High Energy Scattering in 
QCD, hep-ph/0303204. Published in Quark-Gluon Plasma 3, Eds. R. C. Hwa and X.-N. 
Wang, World Scientific, 2003. 

[13] J. W. Cronin, et al. Phys. Rev. Dll (1975) 3105; D. Antreasyan, et al. Phys. Rev. Lett. 38 
(1977) 112; D. Antreasyan, et al. Phys. Rev. D19 (1979) 764. 

[14] Yu.V. Kovchegov and A.H. Mueller, Nucl. Phys. B529 (1998), 451. 

[15] A. Accardi and M. Gyulassy, "Cronin effect vs. geometrical shadowing in d + Au collisions 
at RHIC," arXiv:nucl-th/0308029. 

[16] D. E. Kharzeev, E. Levin, and L. McLerran, Phys. Lett. B561 (2003) 93. 

[17] B.Z. Kopeliovich, J. Nemchik, A. Schafer, and A.V. Tarasov, Phys. Rev. Lett. 88 (2002) 
232303; A. Accardi, "Cronin effect in proton nucleus collisions: A survey of theoretical 
models," arXiv:hep-ph/0212148; I. Vitev, Phys. Lett. B562 (2003) 36; X. N. Wang, Phys. 
Lett. B565 (2003) 116; R. C. Hwa and C. B. Yang, "Final-state interaction as the origin of 
the Cronin effect, " arXiv:nucl-th/0403001. 

[18] For a review, see M. Gyulassy, I. Vitev, X. N. Wang and B. W. Zhang, "Jet quenching and 
radiative energy loss in dense nuclear matter," arXiv:nucl-th/0302077. Published in Quark 
Gluon Plasma 3, editors: R.C. Hwa and X.N. Wang, World Scientific, Singapore. 

[19] M. Braun, Eur. Phys. J. C16 (2000) 337; Phys. Lett. B483 (2000) 105. 

[20] Yu. V. Kovchegov and K. Tuchin, Phys. Rev. D65 (2002) 074026. 

[21] A. Dumitru and L. McLerran, Nucl. Phys. A700 (2002) 492. 

[22] J. P. Blaizot, F. Gelis, R. Venugopalan, "High energy p A collisions in the color glass 
condensate approach. I: Gluon production and the Cronin effect,", hep-ph/0402256; "High 
energy p A collisions in the color glass condensate approach. II: Quark production, " hep- 
ph/0402257. 

[23] R. Baier, A. Kovner and U. A. Wiedemann, Phys. Rev. D 68 (2003) 054009. 

[24] D. Kharzeev, Yu. V. Kovchegov, and K. Tuchin, Phys. Rev. D66 (2003) 094013. 

[25] J. L. Albacete, N. Armesto, A. Kovner, C. A. Salgado and U. A. Wiedemann, Phys. Rev. 
Lett. 92 (2004) 082001. 



85 



[26] E. Iancu, A. Leonidov and L. McLerran, Nucl. Phys. A692 (2001) 583; Phys. Lett. B510 
(2001) 133; E. Ferreiro, E. Iancu, A. Leonidov and L. McLerran, Nucl. Phys. A703 (2002) 
489. 

[27] D. Kharzeev, E. Levin and M. Nardi, Nucl. Phys. A730 (2004) 448, erratum in arXiv:hep- 
ph/0212316. 

[28] J. Jalilian-Marian, "Forward Rapidity Hadron Production in Deuteron Gold Collisions from 
Valence Quarks," arXiv:nucl-th/0402080. 

[29] Yu.V. Kovchegov, Phys. Rev. D54 (1996), 5463; Phys. Rev. D55 (1997), 5445. 

[30] J. Jalilian-Marian, A. Kovner, L. McLerran, H. Weigert, Phys. Rev. D55 (1997) 5414. 

[31] J. Jalilian-Marian, A. Kovner, A. Leonidov and H. Weigert, Nucl. Phys. B504 (1997) 415; 
Phys. Rev. D59 (1999) 014014. 

[32] I. Balitsky, Nucl. Phys. B463 (1996) 99; Phys. Rev. Lett. 81 (1998) 2024; Phys. Lett. B518 
(2001) 235; High-energy QCD and Wilson lines, hep-ph/0101042. 

[33] Yu. V. Kovchegov, Phys. Rev. D60 (1999), 034008; ibid. D61 (2000) 074018. 

[34] H. Weigert, Nucl. Phys. A703 (2002) 823. 

[35] E. Iancu and L. McLerran, Phys. Lett. B510 (2001) 145. 

[36] E. Iancu, K. Itakura and L. McLerran, Nucl. Phys. A724 (2003) 181. 

[37] A. H. Mueller, Nucl. Phys. B558 (1999) 285. 

[38] E. Iancu, K. Itakura, and L. McLerran, Nucl. Phys. A708 (2002) 327. 
[39] A. H. Mueller and D.N. Triantafyllopoulos, Nucl. Phys. B640 (2002) 331. 
[40] D.N. Triantafyllopoulos, Nucl. Phys. B648 (2003) 293. 

[41] S. Munier and R. Peschanski, Phys. Rev. Lett. 91 (2003) 232001 [arXiv:hep-ph/0309177]; 
Phys. Rev. D69 (2004) 034008 [arXiv:hep-ph/0310357]; "Universality and tree structure of 
high energy QCD," hep-ph/0401215. 

[42] A. H. Mueller, Nucl. Phys. A724 (2003) 223. 

[43] L.N. Lipatov, Sou. J. Nucl. Phys. 23 (1976) 338; E.A. Kuraev, L.N. Lipatov and V.S. Fadin, 
Zh. Eksp. Teor. Fiz 72, 3 (1977) {Sov. Phys. JETP 45 (1977) 199); Ya.Ya. Balitsky and 
L.N. Lipatov, Sov. J. Nucl. Phys. 28 (1978) 822. 

[44] V.S. Fadin and L.N. Lipatov, Phys. Lett. B429 (1998) 127; G. Camici and M. Ciafaloni, 
Phys. Lett. B430 (1998) 349. 

[45] G.P. Salam, JHEP 9807 (1998) 19; M. Ciafaloni, D. Colferai, Phys. Lett. B452 (1999) 
372; M. Ciafaloni, D. Colferai, and G.P. Salam, Phys. Rev. D60 (1999) 114036. 



86 



[46] K. Rummukainen and H. Weigert, "Universal features of JIMWLK and BK evolution at 
small x," arXiv:hep-ph/0309306. 

[47] E. Iancu and A.H. Mueller, Nucl. Phys. A730 (2004) 494. 

[48] A.H. Mueller and A.I. Shoshi, "Small-x physics beyond the Kovchegov equation, " arXiv:hep- 
ph/0402193. 

[49] N. Armesto and M. Braun, Eur. Phys. J. C20 (2001) 517; ibid. C22 (2001) 351. 

[50] K. Golec-Biernat, L. Motyka, and A.M. Stasto, Phys. Rev. D65 (2002) 074037. 

[51] E. Levin and M. Lublinsky, Phys. Lett. B521 (2001) 233; Eur. Phys. J. C22 (2002) 647; 
M. Lublinsky, Eur. Phys. J. C21 (2001) 513. 

[52] K. Golec-Biernat and A.M. Stasto, Nucl. Phys. B668 (2003) 345. 

[53] F. Gelis and J. Jalilian-Marian, Phys. Rev. D67 (2003) 074019. 

[54] J. Jalilian-Marian, Y. Nara and R. Venugopalan, Phys. Lett. B577 (2003) 54. 

[55] E. Levin and K. Tuchin, Nucl. Phys. B573 (2000) 833; Nucl. Phys. A691 (2001) 779; Nucl. 
Phys. A693 (2001) 787. 

[56] A.M. Stasto, K. Golec-Biernat, and J. Kwiecinski, Phys. Rev. Lett. 86 (2001) 596. 

[57] V.N. Gribov and L.N. Lipatov, Sov. Journ. Nucl. Phys. 15 (1972) 438; G. Altarelli and G. 
Parisi, Nucl. Phys.B126 (1977) 298; Yu. L. Dokshitzer, Sov. Phys. JETP 46 (1977) 641. 

[58] A. H. Mueller, Nucl. Phys.B643 (2002) 501. 

[59] K. Golec-Biernat and M. Wiisthoff, Phys. Rev. D59 (1999) 014017; ibid. D60 (1999) 114023. 

[60] E. Iancu, K. Itakura and S. Munier, "Saturation and BFKL dynamics in the HERA data at 
small x," arXiv:hep-ph/0310338, to appear in Phys. Lett. B. 



87 



